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Abstract 

The  theory  of  Massera  and  Schaffer  relating  the  existence 
of  unique  almost  periodic  solutions  of  an  inhomogeneous  linear 
equation  to  an  exponential  dichotomy  for  the  homogeneous 
equation  has  been  completely  extended  to  discretizations  by 
a  strongly  stable  difference  scheme.   In  addition  it  has  been 
shown  that  the  almost  periodic  sequence  solution  will  converge 
to  the  differential  equation  solution  at  a  rate  O(k^)  where  p 
is  the  accuracy  of  the  scheme,  uniformly  in  t,  if  the  coeffi- 
cients are  sufficiently  smooth. 

The  preceding  theory  has  also  been  applied  to  a  class 
of  exponentially  stable  partial  differential  equations  to 
which  one  can  apply  the  Hille-Yoshida  Theorem,   It  is  possible 
to  prove  the  existence  of  unique  almost  periodic  solutions  of 
the  inhomogeneous  equation   which  can  be  approximated  by 
almost  periodic  sequences  which  are  the  solutions  to  appropri- 
ate discretizations.   Two  methods  of  discretizations  are 
discussed;  the  strongly  stable  scheme  described  above  and 
the  Lax-Wendroff  scheme. 
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Introduction 

This  work  extends  certain  facets  of  the  theory  of 
Massera  and  Schaffer  [1]  (in  the  future  referred  to  as  M&S) 
relating  properties  of  the  solution  of  an  inhomogeneous 
ordinary  differential  equation  (ODE)  to  properties  of  the 
solutions  of  the  homogeneous  system.   In  Part  I  we  extend 
this  theory  to  difference  approximations  of  the  ODE.  In 
Part  II  we  use  the  ODE  theory,  exclusively  in  the  stable 
case,  together  with  the  Hille-Yoshida  Theorem  to  obtain 
results  for  partial  differential  equations. 

The  study  of  admissibility  theory  and  dichotomy  theory 
for  difference  equations  was  first  done  by  Coffman  and 
Schaffer  [2].   This  work  differs  from  their  work  in  that  we 
are  concerned  with  the  preservation  of  admissibility  and 
dichotomy  properties  by  a  difference  equation  which  is  used 
to  approximate  an  ODE  system  which  has  certain  of  these 
properties. 

First  let  us  review  the  theory  of  M&S.   Consider  the 
equations 

(0.1)  y  =  A(t)y 

(0.2)  y  =  A(t)y  +  f(t) 

Here  the  independent  variable   t   ranges  over  the  whole  real 
line  and  for  each   t   the  vector   y(t)   lies  in  a  Banach 
space   E   (which  may  be  infinite  dimensional).   A(t),  for 
each  fixed  t,   belongs  to  the  space  of  bounded  operators  on  E 
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which  we  will  denote  by  L(E,E).  Let  B  and  D  be  two 
Banach  spaces  of  a  function  from  the  real  line  into  E. 
We  say  the  pair   (B,D)   is  admissible  for  (0.2)  if 

a)  Vf   G   B     3a   unique   solution      y    (to    (0.2))    e    d,    and 
(0.3) 

b)  Ityllj^   <   Kllfllg    . 

This  definition  is  more  restrictive  than  that  of  M&S .   We  are 
not  going  to  deal  with  all  the  subtleties  of  their  theory, 
but  just  with  certain  important  parts.   It  is  pointed  out, 
however,  that  (0.3b)  is  actually  superfluous  (see  M&S, 
Chapter  5)   but  we  include  it  in  the  definition  of  admissi- 
bility for  simplicity. 

Observe  that  (0.3b)  states  that  if  we  write  the   y 
given  in  (0.3a)  as 

(0.4)  y  =  C(f) 

then  C,  which  is  obviously  a  linear  operator  from  B  into  D, 
is  also  bounded.   We  will  be  concerned  almost  exclusively 
with  the  case  that   D  =  B  =  A(E)   where   A   is  the  space  of 
almost  periodic  (AP)  functions  with  range  in   E.   The  range 
in  this  notation  will  generally  be  omitted  and  we  will  simply 
write  A,  etc.,  if  no  confusion  might  arise.   At  times  we 
will  also  be  concerned  with  the  cases   B  =  D  =  L^   or 
B  =  D  =  C   which  we  define  as  the  space  of  bounded  continu- 
ous  functions  with  range  in   E.   Observe  that   A   and   C^ 
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are  closed  subspaces  of  L   and  we  have  the  inclusions 

■>-^  oo 

(0.5)  A  C  C   C  L 

^         '  OO  OO 

A(t)  will  always  be  an  AP  or  L^  operator  function 
unless  stated  otherwise. 

Associated  with  the  concept  of  admissibility  we  define 
the  concept  of  an  exponential  dichotomy  for  (0.1).   Specifi- 
cally (0.1)  has  an  exponential  dichotomy  if  there  exist 
projections   P,,  P2=  I  -  P,   such  that  if  Y(t)  is  the 
fundamental  operator  solution  to  (0.1)   (Y(0)  =  I,   see  M&S , 
Chapter  3) ,  the  following  estimates  hold  for  some  a  >  0. 


(0.6) 


a)  llY(t)P^Y  -^(3)11  <  Ke""^^  ^^  ,     t<s 

b)  llY(t)P2Y""'"(s)ll  £  Ke~"^^~^^  ,     s>  t 


Note  that  this  is  not  the  general  definition  of  an 
exponential  dichotomy  given  in  M&S  (Chapter  8)   but  it  is 
equivalent  to  their  definition  when   E   is  finite  dimensional. 
In  the  case  that  E   is  infinite  dimensional  we  will  in 
general  only  be  concerned  with  the  stable  case;  that  is, 

P;L  =  I'  ^2    =    °- 

Observe  that   P,   is  merely  the  projection  on   E^  ,  the 

space  of  initial  data  of  solutions  to  (0.1)  which  are  bounded 

for   t  G  [0,").   To  see  this  note  that  if  y^   were  in  the 

range  of   P^   and  if  the  corresponding  solution   y(s)   were 

bounded  for   s  >^  0   then,  by  setting   t  =  0  in  (0.6b),  we 

obtain 
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(0.7)        ^0  "  ^  """(3)7(3)  =  P2Y  "^(3)y(3) 


and  Yq  =  0   follow3  from  letting   s  ->  <».   Similarly  P2  is 
ju3t  the  projection  on   E_  ,  the  space  of  initial  data  which 
are  bounded  for   t  G  (-°°/0].   An  exponential  dichotomy  merely 
states  that  E  =  E   ©  E_   and  that  the  solutions,  in  addition 
to  being  bounded,  decay  exponentially. 

Now  it  is  a  fundamental  result  of  M&S  that  if  A(t)  ^  A 
then   (A, A)   is  admissible  for  (0.2)  iff  (0.1)  has  an  expon- 
ential dichotomy.   The  preceding  statement  is  also  true  if 
A  is  replaced  everywhere  by   L_^.   We  prove  here  the  easy  part 
of  the  result;  namely  if  we  have  an  exponential  dichotomy, 
then  we  have  admissibility.   If   f  e  l   define 


(0.8)        y(t)  =  C(f)  =   I  G(t,s)  f(s)  ds 
where 


(0.9) 


G(t,s)  =   y(t)P^Y~"'"(s)  ,       t  >  s 


=  -Y(t)P2Y  """(s)  ,       s  >  t 


Observe  that  the  estimate 

(0.10)  llG(t,s)ll  <  Ke""'^"^' 

holds  by  the  definition  of  an  exponential  dichotomy. 

Using  (0.10)  we  see  that  the  integral  in  (0.8)  exists, 
that   y   satisfies  (0.2),  and  that   y   is  the  unique   L^ 
solution  to  (0.2)  and  in  fact 
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(0.11)  llyll   <  —  II  fll  . 

So  that  C,  defined  in  (0.8),  is  a  bounded  operator  from 
L  (E)  -^  L  (E)  .  To  see  that   C:  A  ^  A  we  merely  let   T 
be  a  common  e-almost  period  for  f  and  A.  Then  w(t)  -   y(t+T) 
-  y(t)   is  the  unique   L^   solution  to  (0.2)   with 
inhomogeneous  term 

(0.12)   [A(t+T)-A(t)  ]y  (t+T)  +  f(t+T)  -  f(t)  =  0(e) 

We    then   use    (0.11)    to   obtain 

(0.13)  II  (y(t+T)    -   y(t))ll^  =  0(e)    . 

The  proof  of  the  converse  is  more  difficult  and  can  be 
found  in  M&S,  Chapter  10.  A  simplified  proof  valid  only  in 
the  finite  dimensional  case  is  given  in  the  Appendix. 

In  Section  1  we  will  define  almost  periodic  sequences, 
which  will  be  the  type  of  solution  we  will  be  searching  for. 
Section  2  is  the  most  important  of  this  work.   Here,  after 
discussing  the  properties  of  the  strongly  stable  difference 
schemes  we  will  be  using,  we  will  introduce  a  transformation 
which  will  separate  out  the  roots  of  the  scheme  inside  the 
unit  circle  and  permit  us  to  work  with  the  Euler  1-step  scheme. 
This  technique  was  originally  developed  by  Engquist  [3] 
although  the  author  was  not  aware  of  his  work  when  the 
formulation  given  in  Section  2  was  developed.   Engquist' s 
results  will  be  discussed  more  thoroughly   at  the  end  of 
Section  3. 
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In  Section  3  we  settle  the  question  of  admissibility 
when  the  homogeneous  system  is  exponentially  stable.   Using 
an  inequality  which  is  an  exact  discrete  analogue  to  the 
Gronwall  inequality,  we  will  show  that  the  homogeneous 
difference  equation  is  also  exponentially  stable.   From 
there  it  will  be  a  simple  matter  to  obtain,  for  any  suffi- 
ciently small  time  step   k,   the  existence  of  an  almost 
periodic  sequence  as  a  solution  to  the  inhomogeneous 
difference  equation.   Furthermore  the  sequences  converge 
uniformly  to  the  unique  AP  solution  to  the  ODE  with  a 
uniform  error  O(k^)   where   p   is  the  order  of  accuracy  of 
the  scheme.   We  will  also  show  that  this  solution  can  in 
fact  be  calculated,  i.e.   it  is  stable  under   roundoff 
errors  and  errors  in  initial  data. 

In  Sections  4  and  5   we  deal  with  the  case  that  the 
homogeneous  system  has  a  general  exponential  dichotomy. 
We  will  show  that  the  corresponding  inhomogeneous  difference 
equation  also  has  an  exponential  dichotomy.   This  is  only 
of  theoretical  interest  as  the  solution  will  no  longer  be 
stable  under  roundoff  errors  or  errors  in  initial  data. 
In  Section  6  we  will  deal  with  some  miscellaneous  topics, 
especially  the  convergence  of  the  mean  value  of  the  AP 
sequence  to  the  mean  value  of  the  AP  solution,  and  also  the 
weakly   nonlinear  case. 

In  Part  II  we  extend  this  theory  to  a  simple  class  of 
partial  differential  equations  which  can  be  written  as  an 
evolution  equation 
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(0.14)  y  =  [B-6]y  +  f 

where   B   is  an  unbounded  operator  which  satisfies  the 
conditions  of  the  Hille-Yoshida  Theorem  and   6 (t)   is  an 
AP  function  such  that  the  homogeneous  system  is  exponentially 
stable.   We  can  obtain  a  unique  AP  solution  to  (0.14)  by 
using  the  formula  which  would  be  valid  if  B  were  bounded 
and  then  showing  that  under  mild  restrictions  on   f   the 
resultant  function  does  in  fact  satisfy  (0.14). 

In  Section  8  we  introduce  a  family  of  bounded  operators 
B,   which  are  spatial  discretizations  to   B.   We  construct 
functions   y,   which  are  the  unique  AP  solutions  to 

(0.15)  yj^  =  [Bj^-6]y^  +  f  . 

We   will    give   conditions    to   insure    that 

(0.16)  II  (y-yi^)ll^=  O(h^)    , 

where   j   is  the  order  of  the  approximation  of  B,  to  B. 
In  Section  9  we  apply  the  theory  of  Part  I  to  (0.15)  to 
obtain  an  AP  sequence   y   ,   which  approximates   y,  . 
We  have  however  the  unfortunate  restriction 

(0.17)  rk=o(^) 

h 
where      IIB.II    =   0(l/h    ).      A  more    favorable    result 

(0.18)  ^  =   0(1) 

h 

is  obtained  in  Section  10   when  using  the  Lax-Wendroff  scheme 
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In  Section  11  we  will  consider  the  extension  of  these 
results  in  the  case  that  the  operator  B   is  perturbed  by 
some  bounded  AP  perturbation  D(t). 

Finally  in  the  Appendix  we  will  give  a  proof  of  the 
basic  ODE  theorem  that  admissibility  is  equivalent  to  the 
existence  of  an  exponential  dichotomy.   The  proof  is  valid 
only  in  the  finite  dimensional  case  but  is  simpler  than  the 
proof  given  in  M&S  and  is  also  simpler  than  a  finite 
dimensional  proof  to  be  found  in  Coppel  [4], 
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1 .   Almost  Periodic  Sequences 

Our  first  task  is  to  introduce  the  discrete  analogue 
of  an  AP  function.   Following  Corduneanu  [5]  we  define  an 
almost  periodic  (AP)  sequence   a    with  range  in  E  as 
follows : 

a   is  AP   iff  given   e  >  0  there  exists  a  length  L(£) 
n  ^ 

(a  positive  integer)  such  that  in  any  sequence  of  L 
consecutive  integers  there  exists  an  N  such  that 

(1-1)  "%+N-^n"c»<^ 

where  the  sup  in  (1.1)  is  taken  over  n.   As  shown  in 
Corduneanu  (page  45)  this  is  equivalent  to  normality  i.e. 
given  any  sequence  of  integers   N^   the  sequence 

b   .  =  a   ,  will  have  a  uniformly  convergent  subsequence. 

n ;  1    n+N . 

Although  the  proof  given  in  Corduneanu  is  stated  for  scalar 
valued  sequences  this  proof  is  obviously  valid  if  the  range 
is  any  Banach  space  E. 

If  we  define  the  space   L^(E)   as  the  Banach  space  of 
bounded  sequences  with  range  in  E  then  the  AP  sequences  A  (E) 
form  a  closed  subspace  of  l"(E).   As  usual  the  argument  E 
will  be  omitted  when  no  confusion  can  arise. 

We  point  out  that  normality  can  be  used,  exactly  as  in 
the  continuous  case,  to  show  that  for  any  finite  set  of  AP 
sequences   a''",...,aj^  with  range  in  possibly  different  spaces 
E"'",...,E^   and  for  any   e  >  0,  there  is  always  a  length  L(£) 
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such    that   in   any   interval   of    length   L  we   can    find    a  common 
e-almost  period. 

Finally  we   note    that    if      f (t)    is    an   AP    function   then 
the   sequence      f      =    f(nk)       is    an   AP   sequence    for   any    real   k, 
The    converse    is    also   true    as    shown   by   Corduneanu    (page    47) 
but  we    shall   not   use    that. 
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2 .      Properties    of   the   Difference   Scheme 

We    consider   linear      il-step,    strongly   stable   schemes 
described   as    follows. 

(2.1)  I      a.y       .=k      J      3.y       .=k      f      6.[A^.y_^.  +  f^.] 
j=0      =>    """^^  j=0      :i    "+3  j  =  0    ^       n+3^n+D       n^j 

Here  we  assume  that  we  are  discretizing  the  ODE 

(2.2)  y  =  A(t)y  +  f(t) 

and     y      =   y(nk),      A     =   A(nk),       f     =    f(nk)      where      k      is    the 
-^n        -^  n  n 

time   step. 

Associated  with  (2.1)  we  have  the  polynomials 

i  .  i  . 

(2.3)  p(x)  =  I      a.x^  ,    a(x)  =  I      6  .x^  . 

j=0   ^  j=0   ^ 

It  is  well  known  (see  Dahlquist  [6]  or  Heinrici  [7])  that 
consistency  implies  that   x  =  1   is  a  simple  root  of  p(x)  =0 
and  that 

(2.4)  p'  (1)  =  a(l)  =1 

where  we  have  normalized  the  coefficients  so  that  the  common 
value  in  (2.4)  is  one. 

Strong  stability  of  the  scheme  is  achieved  by  restricting 
the  size  of  the  other  £-1  roots  of  p(x)  =  0.   Specifically, 
if  we  number  these  roots    x   ,  u  =  l,...,l,      then  we  require 
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that   there  be   a  positive  number      0    <    1,    such   that 


(2.5)  X,    =1,  |x|<e<l,  u=   2,...,£    . 


Finally,    we   make    two    further   assumptions, 


(2.6) 


X      distinct    / 
u 


U.     —"     ^  f  •  •  •   f  Xj      f 


(2.7) 


X,      ^      0  , 


U     ~~     ^/•••/X' 


(2.7)  is  necessary  because  to  get  an  AP  solution,  the 
difference  equation  should  be  solved  backwards  and  forwards. 
This  assumption  can  be  removed  in  the  important  case  that 
the  homogeneous  ODE  system  is  exponentially  stable.  Condition 

(2.6)  can  be  removed  in  all  cases  and  is  included  here  only 
to  simplify  the  following  proofs.   The  removal  of  these 
conditions  will  be  discussed  in  Section  6. 

Now  in  working  with  a  multistep   scheme,  the  standard 
procedure  is  to  convert  it  into  a  one-step  scheme.   To  do 
this  we  define  the  space   E„  =  ExEx...xE.   (We  will  usually 
write  vectors  in  E,  in  coliomn  vector  form.)   We  give   E.  the 
norm  inherited  from  this  definition,  namely  if  w  e  E„  and 


£ 


^  y' 


(2.8) 


then 


(2.9) 


w 


wil       =   max   II  y    II 
i 


Here  we  have  explicitly  indicated  the  E  norm  in  (2.9). 
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We  now  consider  the  discretization  of  the  linear 


inhomogeneous  system  (2.2).   If   w  = 


^n+£-l 

^  y' 


then  we  get 


n 


(2.10) 


w  ,  T  =  C  w   +  kf 
n+1    n  n     n 


Here  if  we  define   z   =  k  A   , 

n      n 


(2.11)   C^  =  C(z^ z^^^)  = 


(V-^£^n+£)"'(^il-l^n+£-r'^£-l^)  '  '  •  ^  V'^il^n-H^)  "' ^^0  V«0^^ 


I 
0 


0 

I 


f  = 

n 


(a,I 


Observe  that   C   ^  L(E„,EJ  and  will  be  just  an  m£xm£ 
n       £   £ 

matrix  in  the  case  that  E  is  an  m-dimensional  space  (in  which 
case  it  is  called  the  companion  matrix) .   Note  also  that  since 
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A   is    uniformly  bounded   in      t,       the    implicit   term    (a.I-B    z        ) 
can   be   inverted,    for      k      sufficiently   small,    uniformly   in   n. 
We    finally   point   out    that   if  A   and    f  are   AP    (L    )       then   C 
and    f^    are      AP    (L^)       and   also   that      C   given   in    (2.11)    is    a 
smooth    function    (in   the   Frechet   sense)    of   its    l+l   arguments. 
Now   with   A   a   constant   and    z  =   kA  e   L(E,E)      we    consider 
for   small      II  zll       the   homogeneous   difference    scheme 


(2.12)      w^^^   =   U(z)w^ 

U(z)    =    (aj^I-3^z)"^(6^_^z-aj^_^I)     ...     (a^I-3^z)  "^  (6^  z-a^I) 


I 
0 


0 
0 


U(z)  is  a  mapping,  defined  for  small  II  zil  ,  from  L(E,E)  into 
L(Ej,E„).   Observe  the  following  properties  of  U(z). 


(2.13)  (a)  U(z)  is  a  smooth  (in  the  Frechet  sense)  function  of  z 
(b)  U(0)  has  eigenvalues  exactly  x  with  eigenspaces 


E   ,  where 


u 


il-1 

X    y 

u  .  ^ 


X  y 
u  ■' 


},      y  e  E 


(c)  E   = 


1       E 
u=l 


u 


-14- 


Now  (a)  and  (b)  should  be  clear.   To  see  (c)  ,  note  that  the 


spaces   E    are  closed  and  also  that  any  two  have  only  zero 

^1 

in    common.      Let      w  =    .  and    suppose      w     has    an   expansion 


X 


I 


(2.14) 


£ 
w  =       I 
u=l 


r  i-l 
X       y 
u   .    u 


X  y 
u   u 


^  u 
Then  if  we  define  the  vector  w  to  be 

(2.15)  w    =:   V  w 


V   =    X 


l-l 


X 


£-1 
I 


^1^ 


we   can   write 


Xn 


i.e.  V  is  just  the  Vandermonde  matrix  associated  with  the 
distinct  numbers  x,,...,x  .  (2.13c)  now  follows  directly 
from   the    invertibility   of     V. 

The    representations    (2.14)    and    (2.15)    enable    us    to 
define   an   equivalent  norm  on    the   space   E^.       Specifically   if 
w   is    expressed   as    in    (2.14)    we    define 


(2.16) 


wi 


u 


max   II  y    II  „ 
-^  u  E 


u 


n  £ 

That   II  II   is  equivalent  to   II  11    follows  immediately  from 

the  representation  (2.15)  and  the  invertibility  of  V.   Note 

chat  the  equivalence  of  these  norms  implies  the  equivalence 

of  the  operator  norms  they  induce  on   L(E^,E^).   Of  course 
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this  paragraph  is  superfluous  in  the  case  that  E  is  finite 
dimensional. 

For  future  use  we  point  out  that  (2.13b,c)  imply- 
that  U   (z)   exists  for  small   II  zll   and  is  smooth  in  z. 

1       r  1 

Now  define  the  space  E   =  ©   )   E  .   We  have  E  =  E,©E 

u>l   ^  ^1 

and  we  note  that  E-,  is  canonically  isomorphic  to  E.   At  times 

we  will  identify  E,  with  E  but  this  should  not  cause  any 

confusion.   With  respect  to  the  decomposition  E  =  E,©  E  we 

see  that   U(0)   is  in  block  diagonal  form;  symbolically, 

rD(0)    0 
(2.17)  U(0)  =1 

^   0   B(0) 

here      D(0):    E,    ^   E-,       and   is    the    identity,    while   B(0):    E   ^E 

and   II B  (0)11       <    0.    This    is    the    operator  norm  induced   on 

L(E„,E„)    by    the    II     II       norm  on   E,      and    follows    from   the    fact 

that   B(0)    is    just   multiplication   by   x      on    the   space   E    . 

We    can  now   state    the    fundamental    theorem  of   this    section. 


Theorem   1.       For   small      II  zll       there   exists    an   operator 
T(z):    L(E,E)     ->    L(E„,E^)     such    that 

(a)  T(0)    =    I       (Identity   on   E.) 

(b)  T(z),    T       (z)       are    smooth   in    z 

(c)  L(z)    =   T       (z)U(z)T(z)     is    in   block    diagonal    form 
with   respect   to   the   decomposition      E.    =   E,    ©   E    . 

(d)  Writing  L(z)    symbolically    as 

■D(z)  0 


L(z)    = 

0         B(z)    ' 
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then  B(z)|   .  =  B(0)  (from  (2.17))  and  llB(z)ll"  <  G  while 
D(z)  has  an  expansion 

(2.18)  D(z)  =  I  +  z  +  0(11  zll^) 


where      E,       is    identified  with   E. 

Theorem   1    is    basically    trivial    and   the    proof    involves 
familiar    arguments.      First  we   note    that   U(z)    is    smoothly 
invertible    for   small    II  zH  .       (The    restriction    "for   small    II  zll  " 
will   not  be    stated   explicitly    in    the    future.) 

Next  we   observe    that    (wl    -   U(z))  exists    for   w   in    a 

small    annulus    around   the   circle       |w|    =9    and   this    holds 
uniformly    in    z.      This    follows    from    (2.13b)    and   the    geometric 
series. 

Now   define    the   projections 

(2.19)       P(z)    -    27i      f         (wI-U(z))~"'"    dw    ,         Q(z)    =    I    -   P(z)     . 

I  w|=e 

That  P  and  Q  are  projections  is  a  familiar  result  which 
follows  from  the  resolvent  identity, 

(Wj^I-U(z)  )""""-  (w2l-U(z)  )"■" 

(2.20) 

=  (W2-Wj^)  (w^I-U(z)  )~-'-(w2l-  U(z))~   , 

2 

calculating  P   by  integrating  around  two  slightly  different 

circles,  and  interchanging  the  order  of  integration.  It 
follows  from  the  construction  that  P(z)  and  Q(z)  commute 
with  U(z)  and  that  they  are  smooth  functions  of  z. 
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It  is  clear  by  applying  P(0)   to  an  arbitrary  vector 
w  G  Ej,   and  using  the  expansion  given  in  (2.15)  that  P(0) 
is  exactly  the  projection  onto  E   along  E,  and  Q(0)  is  the 
projection  onto  E,  along  E  . 

Now  define 


(2.21)  T(z)  =  P(z)P(0)  +  Q(z)Q(0)  . 

Clearly   T(0)    =   I    and  T(z)    is    smooth  whence   T       (z)    exists 
and   is    also   smooth.      Note    that    this   holds    in   either   of    the 
norms      II     11^      or      II     11^      on   L(E^,E^).       Let 

(2.22)  L(z)  =  T~'''(z)U(z)T(z)  . 


We  claim  L  is  in  block  diagonal  form  with  respect  to 
the  decomposition   E.  =  E^  ©  E   .   This  follows  directly 
from  the  fact  that   P(z)  and  Q(z)  commute  with  U(z).   In 
fact   let  X  s  E,  ;  then 

(2.23)  U(z)T(z)x  =  U(z)Q(z)x  =  Q(z)U(z)x  =  w   (say). 

Now  if   y  =  T   (z)w   and   y  =  y-,  +  y    with   y,  s  e, 
and  y   G  E   we  have 

(2.24)  w  =  Q(z)U(z)x   =  T(z)y  =  Q(z)y^  +  P(z)y"'- 
whence 

(2.25)  P(z)y"''  =  T(z)y-'-  =  0 

and  so   y   =0   by  the  invertibility  of  T(z).   Thus  L(z) 
maps   E,  "^  E,   and  a  similar  argument  shows  L(z)  :  E   -»■  E  . 
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If  we  write  (using  the  notation  of  (2.17))   L(z)   as 


D(z)     0      D(z)  :  E,  ->  E 
(2.26)     L(z)  =  ^  1: 

0    B(z)     B(z)  :  E   ^  E 


then   we    see    that   D(z)     and   B(z)     are    smooth    in    z;    D(0)    =    I 
and   II B  (0)11       <    9.       It    follows    immediately  by    continuity 
that 


(2.2  7)  llB(z)ll  ^    < 


for    small    II  zll     .      Now    (2.27)    will    imply    that    given    Zw.../Z 


N 


with    II  z  . II    small      then 


^  £         „    .N 


(2.2  8)  II      I     I     B(z  .)  II        <    K 

i=l  ^ 

for   a   constant      K      independent   of   the      z.      and     N. 

It   only    remains    to   study    the    first   order   structure    of 

D(z)    in   order   to   obtain    (2.18).      Expanding   D(z)    about   z   =    0 

we   can  write 

(2.29)  D(z)    =    I    +    C(z)    +    0(11  zll  ^) 

where      C    is    the   Frechet   derivative   of   D(z)    at    z   =    0    and   is 
a  bounded    linear  map    from  L(E,E)    into   L(E, ,E,)    which  we 
identify   with    L(E,E).       It    is    necessary    to   show    that   C    is 
the    identity,    and   as    one   might   expect   this   will    follow   very 
easily    from  consistency. 

Let  A  be   an    arbitrary   element   in   L(E,E)    and   let   z   =   kA, 
a'    =   C (A) .      For      y  e   E      and   identifying  E,    with   E,    we 


-19- 


calculate  the  following  limit 

(2.30)  lim  D(kA)"y  =  (I  +  kA '  +0(k^))"y 
k->-0 

nk=  t 

for  any  fixed   t  >  0. 

tA* 
As  one  would  expect  this  limit  is   e    y   because 

the  difference  scheme  expressed  in  (2.30)  is  consistent 

with  the  ODE   y  =  A'y. 

More  precisely  since 

(2.31)  llD(kA)ll  £  1  +  kR 

for  small  k,  where   R  is  some  fixed  constant,  it  follows 
that 

(2.32)  llD(kA)'^II  <  e^"^  ,     Nk  <  T,   0  <  t  <  T 


Let   x   solve  the  difference  equation 

x  ,  ,  =  D  (kA)  X 

(2.33)  """^         " 

Xq  =  y 

ri  ■"  A.   kn  ^ 

i.e.    X      =   D(kA)    y.      Now   if      x      =   e  y      then      x„      satisfies 

n  -^  n  -^  n 

the  perturbed  difference  equation 

(2.34)  X  ^,  =  [D(kA)  +  O(k^)  ]x^ 

n+1  n 

2 

where  the  term  0(k  )  is  uniform  for  nk  <_   T.  Equation  (2.34) 

follows  directly  from  the  fact  that   x(t)  satisfies  the  ODE 
(2 .  35)  X  =  A  X  . 
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Letting      w      =    x      -    x        we   can    show   easily   enough    that 
^        n  n  n  j  n 


(2.36) 


Iw    II    =    0(k)     , 
n 


nk    <   T    . 


In  fact  from  (2.33)  and  (2.34)  we  obtain 


n 


(2.37) 


w 


=  k    y   D(kA)"  ^  0(l)w.  ^ 
j=l  ^-^ 


If  we    use    (2.32)     to  bound    the    powers    of   D(kA)     and   thus    H  x    II  , 
we   then   obtain    (2.36)    and    letting     k   ->■    0      we   obtain 


(2.38) 


lim   D(kA)    y   =    e         y 
nk=t 


Recalling  the  identification  of  E,  with  E  we  define  the 

_   y 

vector  y  =   •   whence   y  s  e,. 

y 

Thus  (2.38)  can  be  expressed 


r  A  t    -) 
e    y 


(2.39) 


lim  L  (kA)  y  = 
k^O 
nk=  t 


A't 
e    y 


Now  we  are  going  to  show  that  A'  =  A  by  calculating  this 
limit  directly  from  the  definition  of  L (kA) .   In  fact  since 
llu(0)ll^  =  1   it  follows  that   llu(kA)"ll"   is  uniformly  bounded 
for  nk  <  T   (and  of  course  this  is  also  true  if  we  replace 


11    11^  by 


I 


) .      Then  writing 


(2.40)  L(kA)"    y      =      T    """  (kA)     U  (kA) '^    T  (kA)     y    , 


and  letting   k  ^  0   such  that   nk  =  t ,   using  the  continuity 
of   T(z)   and  the  fact  that  T(0)  =  I   together  with  the 
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uniform  boundedness  of  the  powers   U(kA)  ,   we  are  left  with 


At 
e   y 


(2.41)      lim  L{kA)"  y   =   lim  U(kA)'^  y   = 


k^O  k^O  ^^ 

nk  =t  nk=t  e   y 


the  last  equality  follows  from  the  fact  that  a  consistent 
and  stable  difference  scheme  is  convergent. 

Since  A  was  arbitrary,  (2.41)  together  with  (2.39)  shows 
that  A  =  a',   which  proves  that  the  map  C  is  the  identity 
and  this  completes  the  proof  of  Theorem  1. 
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3.   The  Exponentially  Stable  Case 

We  now  assume  that  the  homogeneous  system  (0.1)  is 
exponentially  stable.   Thus,  if  Y(t)  is  the  fundamental 
operator  solution  to  (0.1)   (Y(0)  -    I)  ,   then  the  estimate 


(3.1)  llY(t)Y~"'-(s)ll  <  Ke  "^^  ^^  ,      t  >  s  , 


will  hold  with  positive  constants   K  and  a.   In  the  sequel 
we  will  use   K   as  a  generic  positive  constant  so  it  will 
appear  in  contexts  other  than  (3.1). 

We  wish  to  prove  a  similar  estimate  for  the  homogeneous 
difference  equation  (see  (2.10),  (2.11)) 


(3.2)  w  ^,  =  C  w   . 

n+1     n  n 

Now  C    is  invertible  for  all   n   by  (2.7).   Then  (3.2) 
n 

has  a  fundamental  solution   W    (W.  =  I)   which  is  simply 

n-1 


(3.3) 


W   =  \ \    C.  ,  n  >  0  , 

n  1 

i=0 


-n     _-, 

W    =  ]~T  C.  ,  n  >  0  . 

-n    I  '    1 
i=-l 


It  can  easily  be  verified  that  the  unique   solution 
to  (3.2) ,   w   ,  given  initial  data  w.  is  simply 

(3.4)  w   =  W  w'-'-w.  . 

n     n  D   3 

Note  that  for  n  >  j,   WW.    is  simply 
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-1      ^-1 

(3.5)  W  W.   =  ]  r  C.  ,  n  >  j  . 

In  the  stable  case  we  will  only  use  W  W.    for  n  >  i 

-^  n    ]  —   -■ 

and  we    can   see    that    the   invertibility   of      C-       is   not   required 
in    this    case   and   assumption    (2.1)    is    therefore   not   required 
here. 

We    intend   to   show   that   for   small   k       (this   will   not  be 
stated  explicitly   in    the    future)       (3.1)    implies 

,  -a,k(n-j) 

(3.6)  IIWW.Il£Ke  ,  n>_j, 

where   K   is  used  as  a  generic  constant  and  can  be  taken 

independently  of  k,   while   a,  =  a  +  0(k).   Thus  since  a  >  0, 

a,  can  also  be  taken  independent  of  k   for  small  k. 

To  attain  (3.6)  we  suppose   A(t)   is   C  .   All  derivatives 

that  are  assumed  will  always  be  required  to  be  AP  or  L  .  This 

will  not  be  stated  explicitly  in  the  future.   The  smoothness 

condition  on   A  will  be  removed  in  Section  6. 

Expanding     A    ,  .      about     A        we    can  write    (3.2)    as 
^  ^        n+j  n 

(3.7)  w^,,    =    U(kA„)w^    +    k^O(l)w       . 

n+i  n      n  n 

Here  the  0(1)  term  is  uniform  in  n   and  the  operator   U   is 

given  by  (2.12).   We  are  going  to  show  that  perturbations  of 

2 

the  type   k  0(1)   preserve  exponential  stability.   We  thus 

consider  the  unperturbed  system 


(3.8)  w^^^  =  U(kA^)W^ 
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If  we  let   V  =  T(kA  )w   ,  where  the  operator  T  is 
n        n   n  '^ 

introduced  from  Theorem  1,  and  note  that  by  the  smoothness 

2 

of  T  and  A  we  have  T(kA  ,,)  =  T (kA  )  +  k  0(1),  we  see  that 

n+1        n 

we  should  consider  the  difference  equation 


(3.9)  V  ^,  =  L(kA  )v 

n+1       n   n 

where  the  operator   L  has  been  studied  in  Theorem  1. 

If  V   denotes  the  fundamental  solution  to  (3.9)  then 
n 

since   L   is  block  diagonal  we  see  that  V    itself  is  also 

^  n 

block  diagonal.   If  we  write   V   in  the  form 


(3.10)  V^  = 


r  V^    0 

n 

0   V^ 
n 


we  can  see  that  the  second  part,  the  contribution  from  the 
roots  inside  the  unit  circle,  causes  no  difficulty.   In  fact 
for  n  >  j   we  have 

-1   n-1 


(3.11) 


V^V^    -  I  I  B(kA,) 
n  n      11      1 
-•      1=3 


2..T  1„    ..  ,n-j 


V  V?   II  <  K 
n  3 


where  we  have  used  (2.2  8)  and  K  as  a  generic  constant. 

Since   9   is  a  fixed  number  less  than  1   we  see  immediately 

2 

that  (3.6)  holds  for  V   where  a,   can  in  fact  be  taken 

n  1 

equal  to   a   for  small  k. 

We  must  now  deal  with  V  .  Neglecting  perturbations 

2 
of   k  0(1)   and  identifying  the  space   E^  with   E  we 
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consider  the  difference  equation 

(3.12)  x^^j^  =  (I  +  kAn)x^ 

i.e.      we   have    reduced   the  problem   to    the   study   of   the   Euler 

1-step   scheme. 

We  will   compare      X      ,    the    fundamental   solution    to    (3.12), 

n 

with   Y    (=  Y(nk))   where   Y   is  the  fundamental  solution  to 

n 

the  homogeneous  equation  (0.1). 

Let   H(t,s)  =  Y(t)Y~  (s).   As  a  function  of   t, 
H  satisfies  the  equations 

H  =  A(t)H  , 

(3.13)  H  =  A^(t)H  +  A(t)H, 

H(t,t)  =  I 

Integrating  (3.13)  we  obtain 

1 
H((n+l)k,nk)  =  I  +  k   H(nk+0k,nk)  d0 

0 


(3.14)  ^    ^ 

2 

=    I  +  kA(nk)  +  k 


e  de  dcj)  H(nk  +  (J)ek,  nk) 


0  0 
Now  since   A   is  bounded  it  is  a  standard  result  (M&S, 
Theorem  3.1  C) ,   that 

(3.15)       llH(t,s)ll  <  K^(R)   for   |  t-s  |  £R  , 

where   R   is  any  positive  number.   This  does  not  require 
exponential  stability  and  in  fact  is  an  immediate  consequence 
of  Gronwall's  inequality.   It  now  follows  from  the  second 


-26- 


2 
equation  in  (3.13)  that  the  coefficient  of   k    in  (3.14)  is 

bounded  uniformly  in  n   and  k.   This  implies  that   Y    is 

the  fundamental  solution  operator  to  an  equation 


(3.16)   Y  ^,  =  [I+kD  (k)]Y    =   [I+kA   +  k^O(l)]Y 
n+1         n      n  n  n 


where  the   0(1)   term  is  uniform  in  n.   This  means  that  we 
can  regard  the   x   equation  (3,12)  as  a  perturbation  of  the 
Y    equation.   We  point  out  that  the  reduction  we  have  obtained 
does  not  use  exponential  stability  and  so  is  valid  in  the  case 
that  (0.1)  has  a  general  exponential  dichotomy.   If  we  use 
the  stability,  however,  we  see  that   Y    satisfies  the 
estimate 

(3.17)  IIy^y'^II  <  K  e"''^^''"^^  ,      n  1  j  . 

We  must  now  prove  the  proposition  that  exponential  stability 

2 
is  preserved  under  k   perturbations. 

Theorem  2.      Consider  two  difference  equations 

(3.18)  y^^^  =  R(n,k)y^ 

(3.19)  x^^-,^  =  R(n,k)x^  +  kS(n,k)x^ 

defined   for   small   k.      Suppose    that    (3.18)    is   exponentially 
stable,    that   is      there   exists    constants      K      and      a      independent 
of      k      such    that 

(3.20)  IIY   yT^II     <    K   e'"^^"^"^^     ,  n    >    j    , 

n    j        —  ~ 
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where   Y    is  the  fundamental  solution  to  (3.18).   We  point 

out  that  we  do  not  require   R(n,k)   to  be  invertible  since 

we  only  study  the  solution   Y  Y.    for  n  >^  j .   In  that  case 

Y   is  not  invertible  and  we  must  use  (3.5)  in  place  of  Y  Y .   , 
n  "^         n  D 

but  we  will  retain  this  notation  for  simplicity. 

Under  this  hypothesis  there  exists  an   e_   such  that  if 
for  small  k. 


(3.21)  IIS(n,k)ll^  1  e  1  ^0 

where  the  sup  in  (3.21)  is  over  n,  then  (3.19)  is  exponentially 

stable  and  if   X    is  the  fundamental  solution  to  (3.19)  (the 

n 

dependence  on  k  has  been  suppressed)   the  estimate 

^  ^_i        -a,k(n-j) 

(3.22)  IIX  X.  II  <  K,e 

n  j    —  1 

where  K^  can  be  taken  independent  of  k  and  a,  =  a  +  0(e). 
Before  proving  this  we  note  that  the  case  where  S(n,k) 
=  0(k)  is  automatically  covered.  The  more  general  formula- 
tion will  be  used  in  Section  6  and  more  importantly  in  Part  II. 
Note  also  that  it  is  a  discrete  analogue  to  the  Gronwall 
inequality  and  the  proof  is  in  fact  immediately  suggested  by 
the  proof  of  the  Gronwall  inequality.   This  theorem  is  equivalent 

to  a  lemma  of  Engquist  itself  based  on  a  theorem  of  Strang  [8], 
but  the  proof  given  here  is  simpler  than  Strang's  proof  and 
much  more  suggestive  of  the  Gronwall  inequality. 

To  prove  Theorem  2  we  first  note  that  for  n  >  j  we  have 

(3.23)  X  ■-•-I-  ■-■   ■-■-1  ■  -    "    ■■  ■'-'    -  ~-         "-^ 


X.  =  Y  Y  >  +  k    I         Y  Y  .   S-  T  X,  ,X.- 
n  J     n    3  £=j  +  l   "  ^    ^~    ^~-^    ^ 


■■J 
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This    analogue    of   the   variation   of   constants    formula   can  be 
verified   immediately.       Using    (3.20)    we   get 


(3.24)       IIX   X."-""!! 
n    3 

<   K   e-"'^^"-^^    +   kK   e^^e   "j^   ^-ak(n-£)    i,^   --1,,     ^ 

Here      K      is    the    constant   of    (3.20).       If  we    redefine    K    as 

K  e^'^       (for   small   k)       and  define      v       .    =   e""^"    "^r.^^""""       ^^ 

n  /  J  n  J 


can  write 


n-1 


(3.25)     v^^.  <  K  +  kKe   I   v^  ^  ^  =-      s^_^^.     ,  n  >  j 

Ji— J 

From  (3.25)  we  obtain 


whence 

(3.27)   s   .   <   (1+cKk)"-^  s.  ,   <   K,  e^^'^^"-^)  , 
n,3   -  D  'D   -    1 

where   K^   is  defined  so  that   s.  -IK    independent  of  k 
i  D  '  J     -'■ 

and  j   (from  (3.25)  specialized  to  n  =  j+1  and  the  fact  that 
V.  .  =  1)  . 

Inequality  (3.27)   together  with  (3.25)  and  the  defini- 
tion of   V   .   yields  (3.22)  immediately  and  thus  completes 
the  proof  of  Theorem  2. 

Returning  to  our  specific  case,  it  should  be  clear  that 
after  several  applications  of  Theorem  2  we  can  prove 
exponential  stability  for   v^   defined  as   T(kAn)w^  (see  (3.8) 
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and  the  following  paragraph) •   Since  in  terms  of  fundamental 
solution  operators  we  have 

_  _-l    _i     _  _-l 

(3.28)         WW.   =  T   (kAn)V  V.  T  (kA . ) 
n  J  n  3      ] 

we  obtain  exponential  stability  for  (3.8)  and  another  applica- 
tion of  Theorem  2  yields  exponential  stability  for  the  full 
homogeneous  system  (3.2).   It  is  a  simple  matter  to  go  from 
exponential  stability  to  admissibility  for  the  inhomogeneous 
difference  equation 


(3.29)  w  ^,  =  C  w  +  kg 

n+1    n  n    ^n 

We  first  define  admissibility  in  the  obvious  way. 

The  pair  (L  ,L  )  will  be  admissible  for  (3.29),  for  small 

^  OO  '    OO  '  \  t     f 

k,    iff    for   any   sequence      g      in      L^       (3.29)    has    a   unique   solu- 
tion     w      in      L^      and   this    assignment   is    a  bounded   mapping, 
i.e.    there   exists    a  K    (independent   of  k)    such    that 

(3.30)  llw   II       <   Kll  g   II 

n   °°   —        ^n   <» 

An  entirely  analogous  definition  holds  for  the  admissibility 
of  the  pair  (A  ,A  ) .   As  in  the  ODE  case  we  will  show  (L^,L^) 
is  admissible  and  then  show  that  the  solution  is  AP  if  the 
coefficients  are. 

To  show  (L  ,L  )  admissibility  we  merely  write  the  solution 

n      _i 

(3.  31)  w   =  k   y    W  W.   g.  ,  . 

n      >     n  J    D"! 

-i  =  — oo       •>  -* 

Exponential  stability  implies  the  convergence  of  this  series. 
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which  can  be  immediately  verified  to  be  a  solution  of  (3.29). 
To  show  w  e  L^  we  estimate,  using  (3.6) 

oo 

(3.32)       llw    II        <    llg    II     Kk      J      e""^'^    =    II  g    II       K p   <    II  g    II     K 

n   °°   —      ^n   «>         .^rt  n   °°         ,         -ak   —     ^n   <» 

3=0  1-  e 

(using   x/(l-e      )    ^    1   as    x  -^-   0)    and   this    yields    (3.30)    after 
a   redefinition   of   K. 

To   show   uniqueness    is    trivial,       if     w        were    a  bounded 

^  n 

solution  to  the  homogeneous  equation  (3.2)  we  would  have 
for  n  >^  j , 

(3.  33)  w   =  W  wT"*-  w  . 

n     n  :i    ] 

and     w      =    0      follows    immediately   upon    letting      j    ->■   -0°. 

It   remains    to   show    that       (A    ,A    )       is    admissible.       In 

fact   if      N      is    a    common    e-almost   period      for   C         and      g      , 

^  n  ^n 

we   would  have      w    ,.  -   w        to  be    the    unique      L        solution 

n+N    n  ~i     oo 

to  an  inhomogeneous  equation  with  inhomogeneous  term 

(3.34)         ^  [C  ^-   C  ]w  ^^,  +  g  ^^,  -  g 
k    n+N    n   n+N    ^n+N     n 

and    for    fixed      k    >    0      the    almost  periodicity   of     w        is    a 

^  n 

consequence  of  (3.30). 

We  can  remove  the  factor  1/k  in  the  denominator  by  using 
the  expression  (2.11)  for  C   , 

(3.  35)  C  =  C(kA  ,  .  .  .  ,kA  ^,  )  , 

n       n       n+k 

where   C   is  a  smooth  function  of  its  arguments.  Expanding 

C  ,^-C    using  (3.35)   we  see  that  if   N   is  a  common  e-almost 
n+N   n       -3  V     ' 
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period  for  A   and  f  ,  then  N   will  be  an  0(£:) -almost  period 
'^  n      n  ^ 

of  w   ,  with  the  bound  in  the  0(e)  term  independent  of  k. 

Note  that  these  two  paragraphs  follow  only  from  (L_^,L^) 

admissibility  of  (3.29)  and  do  not  depend  on  the  stability 

of  (0.1). 

The  converqence  of  w   to  the  solution  y   also  follows 
^         n  n 

very  easily.   If  y   is  the  lanique  AP  solution  to  (0.2)  we 
have 

(3.36)     I  o'jYn+j-^  I    ^j^n+j 
D  3 


=  k[  I    a  j 

^      0 


(y  (nk+ekj)  -  y(nk) )  de 


I    B  .  [y(nk+kj)  -  y(nk)  ]] 
J 


=   k  g(k,nk) 


(The   summation   is    from   0   to   il    in   the   above.)       Here    g(k,t)    is 
AP    in   t   for  each   k      and 


(3. 37) 


g(k,t)«       =    o(l)     , 


0    ; 


(3.37)  follows  from  the  uniform  continuity  of  AP  functions  on 


the  whole  real  axis. 


Now  if  w  = 

n 


n+£-l 


then  (3.36)  implies  that  w   is  the 


n 


imique  AP  solution   "  to  the  same  equation  as   w^  except  for 


an  error, 


(3.  38) 


f   g(k,nk) 
0 
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and    (3.37)    together  with    (3.30)    yield 

(3.39)  llw   -w   II       =   o(l)     ,  k   ->   0    . 

n      n   oo 

We    can   obtain    finer   convergence    results   by   imposing   some 
smoothness    conditions    on  A(t)    and   f(t).       In    fact   if    the   scheme 
has    order   of   accuracy     p,    and  A   and    f   have   p+1   derivatives, 
then   y    also  has   p+1   derivatives.       It    is    shown   in  Henrici 
(p.    247)    that 

(3.40)  II    y    a.y  ^.-k   y    B-y  ^-ll    <   k^^"^   Gil  y^^  "^11       , 

3d 

where    the    constant   G   depends    only   on    the   scheme.     (3.40) 
together  with    (3.30)    yield 

(3.41)  llw   -w    II       =   0(kP)     . 

n      n  oo 

We  point  out  that  these  convergence  arguments  are  valid  when- 
ever we  have  admissibility  and  do  not  require  stability. 
Before  leaving  the  stable  case  we  would  like  to  discuss 
problems  relating  to  the  computability  of  the  solution  w  • 

Consider  first  errors  in  initial  data.   Suppose  we 
solve  the  exact  difference  equation  (3.29)  but  use  as  initial 
data  Wq (k)  =  w_(k)  +  e (k)   where  w  (k)   is  the  exact  initial 

data  for  the  AP  sequence  solution  w   and  e (k)  is  bounded  for 

^  n 

small  k,  i.e.   Ile(k)ll  1  e^.   Now  the  solution  w^  ,  which  we 
solve  for,  will  be 

(3.42)  ^n  =  ^n  "*"  ^n®^'^)  ' 

and  we  have 
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(3.43)  llw*-w  II  =  0(e  "^^)  ,        n  >  0  . 

n   n  — 

Now  in  general  we  cannot  expect  e(k)  -»-  0  as  k  -^  0. 
This  would  certainly  occur  if  we  could  obtain  initial  data 
consistent  with  the  initial  data  of  the  solution   y  to  the  ODE. 
However  this  solution  is  known  only  by  an  improper  integral 
involving  not  only  the  given  forcing  function  f(s)  but  also 
the  fundamental  solution  Y(t)  to  the  homogeneous  system 
which  we  could  not  expect  to  know  explicitly  unless  A  were 
a  constant.   Thus  the  most  we  can  assert  is  that,  using  fixed 
initial  data  (say  w-(k)  =0),  if  we  integrate  over  a  sufficiently 
long  interval  we  will  get  the  solution  w   up  to  some  exponen- 
tially decaying  error. 

Now  let  us  consider  round-off  error.   If  we  could 
postulate   that  we  solved  a  perturbed  equation 

(3.44)  w  ^,  =  C  w   +  kf  +  0(k^  ^)     ,        n  >  0,   r  >  0  , 

n+1    n  n      n  — 

where  the  error  is  uniform  in   n   then  we  could  relate  this 

* 

solution  to  w  by 
n  -^ 

(3.45)  w**  =  w*  +  k   y   W  wT"""  O(k^)  ,    n  >  0  , 

n      n      .  ^ ,   n  1 
j  =  l 

and  we   can   estimate    the   error   simply. 


n 


00 


(3.46)       Ilk      I      W   W.-^   0(k^)ll     <    k^Kk      I      e    "^^   =   0  (k^)     . 
j=l      "    ^  r=0 

Let  us  further  examine  the  validity  of  (3.44).  If  we 

include  rounding  errors  in  forming  the  companion  matrix  (or 

operator) ,  together  with  the  error  of  taking  only  a  finite 
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expansion  of  the  term  (a. I  -  S^kA   .)     in  the  implicit 
case  (see  (2.11)),  then  we  could  assume  that  the  linear 
part  of  (3.44)  would  be  replaced  by 

**  r+1 

(3.47)  C    =  C   +  0(k    )  . 

n     n 

r+1 
If  we    also    assume    an   error  0(k        II  w   II  )       in    forming   the 

n 

product  C  w   we  see  we  should  postulate  a  system 
^         n  n  '^ 

**  **  r+1         r+1    ** 

(3.48)w  ^^  =  c  w    +  kg   +  0  (k    )  +  0  (k    « w„  II) 
n+1    n  n      ^n  " 

and   in   order   to   justify    (3.44)    we   must   show   that    the    solution 
to    (3.48)    will   be   bounded    for   n    >_  0 . 
We    can   write    this    solution    as 

**  **  ^  -1 

w      =   " 

"  n     u 

D 

(3.49) 


=   W  w„      +   k      )      W^W .    g .    1 
n  n    0  .^^     n    j      j-l 

+   k[    y      W  w"-'-(0(k^)    +   0(k^llw**    II  )  ] 
.^,      n    1  3-1 

1  =  1  J 


Now  if  we  use  the  fact  that 

(3.50)  k  I      e'"'^''   =   0(1) 

r=0 

and   define      h      =           max  II  w .    II    we   can   estimate      h     by 

"         j  =  0 n        =•  " 

(3.51)  h^   =   0(1)    +   k^  h^   0(1) 

where  the  0(1)  terms  are  independent  of  n.   This  implies, 

** 

for   small   k,    the   boundedness    of      II  w^    II       and   so   justifies    (3.44) 

The    techniques    used   in   Section   2    of  block   diagonalizating 
the    companion      matrix  were    first   developed   by   Engquist 
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although  I  was  not  aware  of  his  work  when  my  formulation  was 
developed. 

Engquist's  procedure  may  appear  to  be  more  complicated 
than  the  procedure  presented  here,  but  this  is  because  he 
proves  a  more  general  result;  namely  that  the  companion 

matrix  can  be  block  diagonalized  to  within  an  error  of 

Pi 
arbitrary  order  (say  0(k   )).   (He  also  does  not  restrict 

himself  to  strongly  stable  schemes,  but  this  imposes  stability 

requirements  on  several  additional  homogeneous  systems.) 

If  one  were  to  consider  his  procedure  restricted  to  a  block 

2 
diagonalization  up  to  0(k  )  then,  while  the  two  formulations 

differ  in  the  lines  of  approach,  they  are  equally  simple. 

Engquist  also  shows  that  the  upper  block  will  agree 
with  the  Taylor  series  expansion  of  (0.1),  up  to  terms  of 
order  p,  if  p  <  p,   is  the  accuracy  of  the  scheme.   The 
proof  of  the  first  order  structure  of  the  upper  block  (the 
term  DCz)  in  Section  2)   is ,  I  believe,  somewhat  simpler 
than  his  proof  restricted  to  the  first  order  term. 

I  would  like  to  point  out  that  in  the  case  we  are 
considering  (an  exponentially  stable  ODE  and  a  strongly  stable 
scheme)   the  extra  fineness  of  the  block  diagonalization  being 
carried  to  order  p+1  will  only  give  the  advantage  of  the 
homogeneous  system  having  a  stability  exponent  -a+0(k") 
instead  of  -a+0(k).   In  particular  in  Engquist's  study  of 
the  uniform  convergence  (for  t  >_  0)  of  the  solutions  of  the 
homogeneous  equation  to  the  solution  to  the  ODE  (see  his 
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Theorem   3,    p.    24-27    and   note    that  what  he    uses    as    a   is 
what  we    call    -a)       one  would    still    obtain    uniform  convergence 
of   O(k^)     (neglecting   round-off   and   errors    in   initial   data) , 
however   the   decay   exponent  would  be   -a+0(k)    instead   of 
-a+0(kP)  . 
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4 .   Admissibility  in  the  Case  of  a 
General  Exponential  Dichotomy 

Here  we  are  concerned  with  admissibility  properties 
in  the  case  that  the  homogeneous  equation  (0.1)  has  a 
general  exponential  dichotomy  as  described  by  (0.6)  with 
P„  7^  0.    The  results  will  be  valid  in  the  infinite  dimen- 
sional case,  although  in  this  case  one  must  bear  in  mind 
that  the  situation  described  by  equations  (0.6)  is  not 
the  most  general  form  of  an  exponential  dichotomy. 

First  of  all  consider  arbitrary  homogeneous  and 
inhomogeneous  difference  equations 

(4.1)  w^^^  =  C^(k)w^ 

(4.2)  w^^^  =  C^(k)w^  +  kg^  . 

Here   C    is  defined  for   k  e  (O/k^]  and  is  AP  in  n  for 
n  u 

fixed  k,  and  is  now  assumed  invertible  so  that  a  fundamental 
solution  operator  W   ,  as  described  by  (3.3),  exists  and 
is  invertible  for  all  n. 

Associated  with  (4.2)  we  have  variation  of  constants 
formulas  in  both  the  forward  and  backward   directions, 

n 


(a)      w   =  W  w„  +  k   y   W  W.-'-g.  T  ,      n  >  0 


(4.5) 


}      =Ww„+k   y   WW.g.^, 
n    n  0      >,   n  ]  j-l 
3-1 

0      _-^ 

/   =Ww„-k    y   WW.g.,, 
n    n  0      •   , 1  n  ]   1-1 
j=n+l    -■   -■ 


(b)      w   =  W  w„  -  k    )   W^W."'"g._,,      n  <  0 


as    can   be   easily   verified. 
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Now  we  will  define  an  exponential  dichotomy  for  (4.1) 
as  follows.  Suppose  there  exist  projections  P, (k) ,  P^ (k) 
=   I    -    P,(k),      such    that    the   estimates 


llW  P,W.^II     <    K   e    "'^^^    ^^     ,  n    >    j 

n    1    j         —  '  —   -• 


(4.6) 


W   P^W-^II     <    K   e-^'^^^-") 
n    2    3        — 


j    1    n    , 


hold  with  K  and  a  independent  of  k.  The  only  difference 

between  this  and  the  ODE  case  is  that  the  projections 

may  depend  on  k. 

If  we  have  an  exponential  dichotomy  then  (L  ,L  )  is 

admissible  for  (4.2).   In  fact  the  unique  bounded  solution 

w   is 
n 

00 

(4.7)  w^(k)  =  k  ^    ^i^i-l 

-i  ^  — OO     J    J 


where 


(4.8) 


3     n  1  ]   '  -  -" 

=  -W  P^wT""",  j  >  n+1 


n  2  j  ' 

Equation  (4.7)  is  of  course  suggested  by  the  ODE  case  and 
is  easily  verified  using  (4.6).   The  admissibility  bound 
(3.30)  can  be  derived  exactly  as  in  the  stable  case  as  can 
(A  ,A  )  admissibility. 

We  next  observe  that  (L  ,L  )  (or  (A  ,A  ))  admissibility 
is  preserved  under  perturbations  of  the  linear  term.  In  fact 
if  we  had  a  system 
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(4.9)  w^,,    =   C    (k)w^    +   kS^(k)w^    +   kg„ 

n+i  n  n  n  n  n 

where,    for   small   k. 


(4.10)  lis    (k)ll       <    e- 

n  oo  —      0 

and      e„      is    some   number   to  be   determined,    then   the    unique   L 

0  '  ^  oo 

solution  can  be  found  by  defining  w   ^  =  0  and  w   . , ,  to  be 

-^  ^        n;0  n;i+l 

the    unique      L^   solution   to 

(4.11)  w    ,T     .^,    =   C   w      .     .    +   kS  w      .    +   kg 

n+l;i+l  n   n,i+l  n   n,i  n 

(for   simplicity    the   dependence   of  C      and   S      on   k   has   been 
suppressed)  .       If      ^n^-i    "^    1'    where   K,    is    the    admissibility 
bound   for    (4.2),    the    contracting  mapping  principle   establishes 
(L    ,L    )    admissibility  with   admissibility  bound 


OO  OO 


^1 
(4.12)  K^    = 


2  (1-^0^1^ 


Of  course  we  will  get  AP  solutions  if  D  +kS   and  g   are  AP . 

^  n    n      ^n 

Now  these  two  principles  certainly  settle  the  question 

of  admissibility  when  the  homogeneous  system  has  a  general 

exponential  dichotomy.   Referring  to  the  reduction  obtained 

in  Section  3  we  see  that  after  a  nonsingular  change  of 

dependent  variable,   v  =  T (kA  )w   ,   v   satisfies  the 
^  n       n   n  '    n 

equation  (see  (3.8)  f f . ) ; 

r  I+D  (k)   0 
(4.13)   v^^.  =  I     " 

n 


V  +  k^  0(l)v  +  kg 
n  n     ^n 


Here   g  =  T (kA  ,^)g.      By  the  boundedness  of  T  and  T   , 
^n       n+1  ^n    ^ 
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admissibility  for  (4.13)  is  equivalent  to  admissibility 
for  (4.2).   But  (4.13)  expresses  the  homogeneous  term  as 
a  perturbation  of  the  system 


(4.14) 


X 


n+1 


I+D    0 
n 

0     B 


X 


n 


n 


with  fundamental  solution 


X   =        n 
n      I  n   T 


1-1 


i=l 


B, 


n  >  0  , 


(4.15) 


f    ^- 


X 


-n 


n 


-n+1 


i  =  0 


bT^ 

1 


n  >  0  , 


and  this  has  an  exponential  dichotomy  with  projections 
independent  of  k   and  in  fact  (after  identifying  E   with  E) 


(4.16) 


P   =  P   ©  P(0) 


P   =  P 
2     2 


Here   P,  and  Py    ^^®  defined  in  (0.6)   while   P(0)   is  the 
projection  onto  E   along  E,  (see  (2.19)  f f ) . 

We  thus  obtain  admissibility  in  this  case  and  we  point 
out  that  the  convergence  proofs  given  in  Section  3  did  not 
require  stability  and  remain  valid  here. 
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5 .   Existence  of  an  Exponential  Dichotomy  for 
the  Homogeneous  Linear  Equation 

We  will  prove  here  that  admissibility  for  the  inhomo- 
geneous  difference  equation 

(5.1)  w  ^,  =  C  w  +  kg 

n+1    n  n     ^n 

implies  an  exponential  dichotomy  for  the  homogeneous  equa- 
tion 

(5.2)  w  , T  =  C  w 

n+1    n  n 

The  proof  is  very  similar  to  a  proof  of  the  analogous  property 
for  the  ODE  case  given  in  the  Appendix.   It  is  a  strictly 
finite  dimensional  proof   and  it  may  then  be  simpler  to  think 
in  terms  of  matrices  rather  than  linear  operators.   We  can 
then  regard  T(z)  as  a  similarity  transformation  which  follows 
an  initial  transformation  putting  U(0)  into  block  diagonal 
form.   The  matrix  L(z)  is  strictly  in  block  diagonal  form 
and  it  is  no  longer  necessary  for  us  to  continually  make  the 
qualification  "identifying  E,  with  E." 


If  we  let   V  -   T (kA  )w   ,  (5.1)  is  transformed  into  (see 
n        n   n 


(4.13)) 


I+kD    0 

(5.3)   v^,,  -  '      "^ 

n+1    ,   .      3 

n 


V  +  k  0(l)v  +  kg 
n  n     ^n 


(For  simplicity  we  will  write  the  inhomogeneous  term  in  (5.13) 

as   g    rather  than   g  . ) 
^n  ^n 
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If  we   call    the    leading   matrix   in    (5.3)       R    (k) ,    then    the 
unperturbed   linear   system 


(5.4)  X  =   R   X 

n+1  n   n 

has  an  exponential  dichotomy  with  projections  P,  and  P„  , 
independent  of  k  (see  (4.15)  and  (4.16)) ,  and  we  must  exhibit 
an  exponential  dichotomy  for  the  homogeneous  version  of  (5.3). 

(5.5)  V  ,  ^  =  R  V   +  k^O(l)  V   . 

n+1    n  n  n 

By  admissibility  (5.5)  can  have  no  L^  solution.  Let 
S  (k)  and  S_(k)  be  the  subspaces  of  initial  data  which  give 
rise  to  solutions  bounded  for  positive  and  negative  n  respec- 
tively.  We  have   S   ^  S_  =  {0}  (suppresssing  the  k  dependence) 
If  we  let   S   be  any  complementary  space  we  have 


(5.6)  E   =  S^  ©  S_  e  S 


with  associated  projections 


(5.7)  I  =  P^  +  P2  +  P3 

Our  first  task  is  to  show   S  =  {0}   i.e.  P^  =  0. 

Define   L    to  be  the  submanifold  of  L    of  sequences 
with  only  finitely  many   nonzero  components.   Let   g   ^  L^ 
and  let   v    be  the  corresponding  unique   L    solution  to 
(5.3).   Then  using  (4.5a)   we  see  that  for  large  positive  n 

we  have 

00 

(^•«)         ^n  =  ^nf^O  ^   '^  .1      ^j'  5j-l^ 
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and  since   v   is  bounded  the  vector  in  the  brackets  lies 
n 

in  S    i.e. 

(a)      P2VQ  =    '   ^   -   "-1 


=  -  k  I      P,  V.   g.  T 
j  =  l   2   ^    ="-1 


(5.9) 

°°      -1 
(b)      P3V0  =  -  k  J      P3  V    g 

-,  =  1 

Similarly  looking  at   v   for  n  ^  -°°   and  using  (4.5b)  we 

obtain 

0      _. 
(a)       P^Vq  =  k    I   P^  V    g 

(5.10) 

0 


(b)       P,v     '    ^   -   "~^ 


Now  if  we  simply  set 


,v_  =  k    y   P^  V~      g.  , 


(5.11) 


^0  =  ^1^3^ 


g^  =  0    ,  i  7^  0  , 


where   z   is  an  arbitrary  vector  in   E.  we  see  immediately 

that   P^  =  0.   We  can  also  see  that  for  any  given   g  g  l 

the  unique   L   solution   v   is  given  by 

00 
(5.12)  ^n  =  ^   ^    ^i  ^i-l 

where   G.   is  given  by  (see  (4.8)) 


,n    .,     .   „-l 
(5.13) 


G.  =   V   P,  V.   ,    n  >  j 
J  n   1   D 

=  -V^  P„  V"-"-  ,    j  ^  n+1 
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It  is  of  course   g'^  which  we  wish  to  prove  is  exponentially 
damped. 

Let  G^      be  the  Green's  function  associated  with  (5,4). 
J 

We  know 

(5.14)  IIg"^!!  <  K  e""^!"""^! 

D   - 

Now   define      h"^   =    g"^   -    g'^    .       If      g .    e   l^      and      v        the   corres- 
3  3  3  :  °°  n 

ponding   solution    to    (5.3),    while      x^      is    the      L^      solution 
to   the    inhomogeneous    version   of    (5.4), 


(5-15)  x^^^   =   R^x^    +    kg^    . 

Then      u     =  v      -    x        is    the   unique   bounded   solution   to 
n  n  n 

(5.16)  u^^^   =    R^u^   +    k^    0(1)    v^    . 


We    thus   have 


(5.17)       u^   =   k      I        h'?    g  ==    k^       I         g"    0(1)    V  . 

-i  =  — oo  J  J  -i=  — 00  -J  -^ 


.2       r         -n 

j  =  -oo  -  -  j: 

If  we    fix   i    and    let      q^_^   =    z,      g^    -    0      for   j    ^   i-1      where 
z   is    an   arbitrary   vector   in   E^^    ,    we   then   obtain,    using    (5.14) 
and   the   admissibility   bound    for    (5.3), 

(5.18)  llH^II     <   K^    ,  for   all   n,    i 

where   K,   is  some  constant  independent  of  k.   (The  restric- 
tion "for  small  k"   is  always  understood.)   The  use  of  (5.14) 
then  implies  the  boundedness  of   G.   i.e. 

(5.19)  llG^il  <  K   , 
where  K^  is  independent  of  k. 
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We  will  now  show  that 


(5.20) 


Pj^(k)  -  P^ll  =  0(k) 


Observe  that  P-,   is  independent  of  k.   To  obtain  (5.20)  we 

note  that   Pt-Pi  =  H-.   If  we  define   g.  ^  =   z,  g.  =  0 
11     0  ^1-1       ^j 

(j  7^  i)  where   z   is  an  arbitrary  vector  we  then  have 
from  (5.12) 


(5.21) 


V.  ,  =  k  G^.  "*"  g.  ,  =  k  G^""""  z  . 
j-1       1   ^1-1      1 


Substituting  into  (5.17)  we  obtain 


(5.22)     u  =  kH"z  =  k^ 
n     1 


I 


J=-C 


g"  0(1)  G^"^ 


Now   using    (5.19)    and    (5.14)    we    see    the   sum   in   brackets    in 
(5.22)    is   0(l/k)    whence   we    conclude    that 


(5.23) 


hJii    =  0(k) 


and  setting  n  =  i  =  0  we  obtain  (5.20). 

Note  that  this  proves  that  the  ranks  of  the  stable  and 
unstable  manifolds  are  unchanged  for  small  k. 

We  can  now  show  quite  easily  that  (5.19)  can  be  replaced 
by  an  exponential  decay  factor. 

For  a  certain  small  positive   e   consider  new  systems 

-ek 
with  the  linear  part  of  (5.3)  and  (5.4)  multiplied  by  e     i.e 


(5.24) 
(5.25) 


v^^T  =  e  ^^[R  +  k^  0(1)  ]v^ 
n+1  n  n 


e      -ek  „    e 
X  ,  ,  =  e     R   X 
n+1         n   n 
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Now  by  Section  4  (see  (4.9)  ff . )   the  inhomogeneous  versions 
of  (5.24)  and  (5.25)  will  both  have  (l",L^)  admissible,  if 
e  is  small  enough.   But  the  new  fundamental  solutions  are 

(5.26)  V^  =  e"""  V^ 

<5.27,  Xj  =  e-^  X„ 

and  clearly  (5.25)  will  have  an  exponential  dichotomy  with 
projections   P-,  and  Py    if  ^  i^  small  enough. 

Now  if   P, (k)  and  P„ (k)  are  the  corresponding  stable 
and  unstable  projections  for  (5.24)  while  G.    is  the 
Green's  function  for  (5,24),  the  same  analysis  that  led  to 
(5.19)  will  yield  for  some  constant   K^  , 


(5.28)  "^j;£"  -  ^3 

In   particular   if      j    >_  n+1,     (5.28)    reduces    to 

(5.29)  IIV     P^   V'.h     <   K,   e"^^^^"""^ 

n      2       3         —      3 

and  one  part  of  the  definition  of  an  exponential  dichotomy 
for  (5.5)  will  have  been  shown  if  we  can  show  Pj^  ^  P]^  ' 
P^  =  P   .  But  this  is  a  simple  consequence  of  (5.20).  In 
fact  if   S^  ,  sf   are  the  new  stable  and  unstable  manifolds 
for  (5.24)  then  (5.26)  yields 

(a)  S   c  s^ 
(5.30) 

(b)  S_  C  S_ 

But    (5.20)     shows    that    the    ranks    of   P^    and   P^^    are    equal 
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(being  equal  to  the  rank  of  P-,)  and  similarly  for  P„  and  P_. 
This  shows  that  the  inclusions  in  (5.30)  are  equalities  and 
that 

(a)  P^  =  P, 

(5.31)  ^  ^ 

(b)  pI  =   P2 

Thus  one  part  of  the  requirements  of  an  exponential  dichotomy 
for  (5.5)  has  been  shown  and  the  other  part  will  follow  on 
replacing  -e  by  +e. 

Lastly  we  point  out  that  an  exponential  dichotomy  for  (5.5) 
implies  one  for  (5.1)  since  in  terms  of  fundamental  solutions 
we  have 

(5.32)  W^  =  T"-'-(kA^)V^  T(kAQ) 

whence  the  projections   P-,{k),  P~(k)  for  (5.2)  are  related  to 
P,  and  Pp  by  a  similarity  transformation 


(5.33) 


(a)    P^(k)  =  T~-^(kAQ)P^(k)T(kAQ) 


(b)    P2(k)  =  T  ■^(kAQ)P2(k)T(kAQ) 
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6.   Miscellaneous  Results 

This  section  is  concerned  with  certain  generalizations 
of  the  theory  developed  in  the  preceding  section.  We  first 
remove  certain  restrictions  which  had  been  imposed  previously 
in  order  to  make  the  exposition  clearer.   In  Section  6A  we 
will  remove  the  restriction  that  A  be  C  while  in  Section  6B 
we  remove  the  restrictions  (2.6)  and  (2.7)  on  the  roots  of 
the  polynomial  p(x)  =  0   which  lie  inside  the  unit  circle. 

We  then  consider  certain  trivial  extensions  of  the 
theory.   In  Section  6C  we  remark  on  the  general  L^  (non-AP) 
case.   In  Section  6D  we  consider  the  convergence  of  the 
mean  value  of  our  sequence  solution  to  the  mean  value  of 
the  solution  of  the  ODE  and  finally  in  Section  6E  we 
consider  the  weakly  nonlinear  case. 
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6A.   Removal  of  the  Differentiability  Condition  on  A 

First  consider  the  case  when  the  homogeneous  equation 

(6.1)  y  =  Ay 
is  exponentially  stable;  i.e. 

(6.2)  llY(t)Y"-'-(s)ll  <  K  e"*^^^"^^  ,     t  >  s, 

where   Y(t)  is  the  fundamental  solution  to  (6.1).   We  assume 

the  A  is  AP  but  is  not  C  .   Now  Theorem  2  does  not  require 

2 
a  perturbation  0 (k  )  but  is  certainly  valid  for  a  perturbation 

k  o(l)  (k  ->•  0)  •   This  leads  us  to  expect  that  the  requirement 

that  A  be  C   can  be  replaced  by  the  uniform  continuity  of  AP 

functions  on  the  whole  real  axis.   This  is  indeed  the  case 

as  one  can  verify  with  little  difficulty. 

In  this  subsection,  we  will  merely  trace  through  the 
proofs  in  Section  3  and  indicate  what  changes  must  be  made 
if  A  is  not  C  . 

Discretizing  (6.1)  we  obtain  the  homogeneous  difference 
equation 

(6.3)  w  ^,  =  C  w 

n+1    n  n 

where  C   is  given  in  (2.11).   Using  the  uniform  continuity  of 
A  we  see  that  (3.7)  can  be  replaced  by 

(6.4)  w^_|_^  =  U(k  A^)w^  +  k  o(l)w^ 

where    the    term   o(l)    is    uniform  in   n    (this   will   not  be   stated 
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explicitly  in  the  future) .  If  we  also  use  the  uniform 
continuity  of  A  to  obtain 

(6.5)         T(k  A  ^,)  =  T(k  A^)  +  k  o(l) 

n+l  n 

we  see  that  Theorem  2  will  yield  exponential  stability  for 

(6.3)  provided  we  can  show  exponential  stability  for  the 
system 


^n+1  =  ( 


I+kA    0 
(6.6)  V  _   -  '      ^ 

■       0      B 
n 


V 

n 


The  lower  block  causes  no  difficulty  and  we  are  left  with 
the  system 

(6.7)  x^^^  =  [I  +  k  A^]x^. 

Finally  we  can  show  that  (6.7)  is  exponentially  stable  by 

using  Theorem  2  to  compare  the  fundamental  solution  X   with 

Y   (=  Y(nk)).   This  requires  some  modification,  as  the  proof 

leading  to  (3.16)  used  the  differentiability  of  A  (see 

(3.13)  f f )  . 

Following  the  notation  of  Section  3,  we  define 

H(t,s)  =  Y(t)Y~  (s)  .   We  can  then  write  (compare  with  (3.14)) 

1 

(6.8)  H((n+l)k,  nk)  =  I  +  kH(nk,nk)  +  k    dO  [H  (nk+ek  ,nk) 

°   -  H(nk,nk)  ] 

=  I  +  kA  +  k  I  de  [A(nk+9k)-A„]  H(nk+6k,nk) 
n  n 


+  k 


1 
f 


d9  A  [H(nk+ek,nk)  -  H(nk,nk)] 
n 
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Finally  we  see  that  if  we  use  (6.2)  together  with  the 
uniform  continuity  of  A  on  the  whole  real  axis  and  the 
equation 

(6.9)  H(t,s)  =  A(t)H(t) 

then  (6.8)  can  be  rewritten  as 


(6.10)         Y^,,  =  [I  +  kA^  +  k  o(l)]  Y„ 
n+1  n  n 

and  Theorem  2  is  immediately  applicable. 

We  thus  have  exponential  stability  for  (6.3)  and  hence 
(A  ,A  )  admissibility  for  the  inhomogeneous  version  of  (6.3), 

Finally  we  point  out  that  the  results  of  Section  4, 
regarding  admissibility  in  the  case  that  (6.1)  has  a  general 
exponential  dichotomy  is  equally  valid  if  A  is  not  C  ,  as 
the  contracting  mapping  principle  (see  (4.9)  ff)  would 
certainly  be  applicable  if  the  perturbation  is  k  o(l).   The 
results  in  Section  5,  however,  use  crucially  the  differenti- 
ability of  A  (in  the  argument  involving  the  deduction  of 
(5.18)  from  (5.17))  and  it  has  not  been  possible  to  extend 
this  result  when  A  is  not  C  . 
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6B.   Removal  of  Restrictions  on  the  Roots   of  p(x)  =  0 

In  Section  2  we  prescribed  two  conditions  on  the  roots 
of  p(x)  =  0  which  lie  inside  the  unit  circle;  namely 


(6.11)  X   distinct  ,      u  =  2,...,£, 

(6.12)  X   7^  0       ,      u  =  2,  .  .  ,£. 

(Recall  that  x,  =  1  and  |x  |  <  9  <  1  for   u  >  1) .   We  will 
now  remove  these  restrictions. 

It  should  be  clear  by  now   that  in  the  exponentially 
stable  case  (6.12)  is  unnecessary.   In  fact  the  fundamental 
solution  to  the  homogeneous  equation  is  used  only  for  n  >_  j  , 
and  the  equation  (3.5)  shows  that  the  invertibility  of  the 
linear  term  C-  is  not  required  for  n  >^  j .   Of  course  the 
notation  WW.    is  no  longer  accurate,  but  except  for  this 
detail   the  results  of  Sections  2  and  3  are  valid  without 
assuming  (6 . 12)  . 

We  deal  next  with  the  restriction  (6.11).  (Note  that 
the  root   x,  =  1  is  always  simple.)   The  removal  of  (6.11) 
in  the  finite  dimensional  case  is  trivial.   In  fact  suppose 
the  root   x   has  multiplicity  q  >  1;  i.e. 

(6.13)  P"^(x^)  =  ...  =  p'^"-^(x^)  =  0 

p'^(x^)  ^  0  . 

It  is  shown  in  Henrici  (p.  214)   that  the  effect  of  (6.13) 
is  that  the  eigenvalue  x  of  U(0)  (see  (2.12))  will  now  have 
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nontrivial  Jordan  blocks.   Thus  U(0)  can  still  be  put  in 
block  diagonal  form  (see  (2.17)) 

I    0 


(6.14) 


U(0)  = 


0   B(0) 

where  B(0)  is  no  longer  diagonalizable ,  but  has  all  of  its 
eigenvalues  bounded  by  0  in  the  norm.   Now  since  B(0)  can 
be  put  in  Jordan  normal  form  with  e  instead  of  1  on  the 
superdiagonal,   for  any   e  >  0,   we  can  certainly  choose 
a  matrix  norm  such  that 


(6.15) 


lie (0)11  <  e 


Since  all  norms  are  equivalent  on  a  finite  dimensional  space, 
the  proof  of  Theorem  1  can  now  proceed  exactly  as  in  Section  2, 

The  case  when  E  is  infinite  dimensional   can  be  handled 
in  exactly  the  same  manner;  however  since  we  can  no  longer 
appeal  to  the  theorems  of  linear  algebra,  we  will  have  to 
carry  out  the  proof  in  more  detail. 

It  is  shown  in  Henrici  (p.  214) ,  that  in  the  scalar  case 
a  basis  for  the  generalized  eigenspace  E   corresponding  to 


the  eigenvalue   x    is  {e   . } ._, 

r        ^ '  3    D"-*- '  •  •  •  » ' 


(6.16) 


r      £-1 
'  x 
r 


'r,l 


where 
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r.D 


r=l 


ft    (^-r)x^-     (^-1) 


ft    (n-r)x(^-(J-l)) 
r=l 


i-1 


fr    (J    -    r)x 
r=l 


J    >    1 


One    can   verify   easily    that   in   fact 


(a)       U(0)    e      ,    =   X  e      , 
r,  1  r    r,  1 


(6.17) 


(b)         U(0)    e       .    =   X     e       .    +   e       •    ^     i 


j    >    !■ 


Also  since  the  vectors   e   .   form  a  new  basis  the  £  x  £ 
matrix  with  columns   e   .  is   nonsingular.   Wow  in  order 
to  have  e  instead  of  1  on  the  superdiagonal,  we  replace  e 
by  e   .  where 


ro 


(6.18) 


j-1  - 
e   .  =  e     e   . 


(see  Bellman  [9],  p.  198).   Here   e  >  0  is  to  be  specified. 
The  matrix  with  columns   e   .   is  still   nonsingular  but 
(6.17b)  will  now  be  replaced  by 


(6.19) 


U (0)  e   .=xe   .+ee   .,. 
r,]     r   r,j       r,3-l 
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Now  the  extension  of  this  to  infinite  dimensional  space 

is  very  simple  and  follows  closely  the  procedure  of  Section  2 

For  any  vector  y  s  e   define  the  vector  e   . (y)  s  e„  by 

r ,  3  ^ 


(6.20)  e^    .(y)    =    e^~^ 


(      3 

n 

r=l 
J 


(£-r)x^"^y 


r=l 


T    (j  +  l-r)x  y 


Equation    (6.19)    clearly   generalizes    to 


(6.21) 


U(0)    e^    .(y)    =   X     e         (y)    +    e    e       .    ,  (y)     ,       j    >    1 
■L/J  ri./j  i,ji 


U(0)    e^^-L(y)    =   x^  e^^^(y) 


We  now  define  the  space   E   .  =  {w  G  e   such  that 

r ,  ]         X, 

w  =  e   . (y)  for  some  y  G  e}.   Clearly  the  space  E   .  is 
closed  for  all  r  and  j.   Furthermore  we  have 


(6.22) 


^£  = 


r,J 


This  follows  by  the  same  argument  as  that  given  in  Section  2 


(see  (2.14)  f f ) .   In  fact  if  w  = 
and  we  had  an  expansion 


X- 


i-  Xc 


is  a  vector  in  E 


a 


(6.23) 


r  x. 


^  H   -> 


-   Z   e   .(y   .: 
^^.   r,3  ^r,] 


r/D 


then  the  {y   .}   would  be  related  to  the  {x.}  by  multiplication 
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1  X 

while   E      is    defined   as      ©        I        E       .    ,    then,    exactly    as 


by   a  nonsingular    5,    x    £   matrix   and    this    establishes     (6.22). 

It   also    follows    from   this    that   if  we    define    a   new  norm  on 

E.    by      II  wll       =   max   II  y       .11       then    II  wll       is   equivalent   to   II  wll 

( max   II  X .  II     )  . 
i  ^   ^ 

Now  if  E,  is  defined  as  E,  -,  (the  root  x,  is  simple) 

I    E   .  , 
r>l,j   ""'^ 
in  Section  2,  we  see  that  U(0)  is  in  block  diagonal  form 

with  respect  to  the  decomposition  E   =  E,  ©  E  ,  and  we  can 

write  (see  (2 . 17) ) 


f   D(0)    0 

(6.24)  U(0)  = 

^   0     B(0) 

where   D(0):  E,  ->  E,  and  is  the  identity,  while  B(0):  E   ^  E 
and  for  e  sufficiently  small  we  have  from  (6.19) 

(6.25)  llB(0)ll"  <  e 


The  proof  of  Theorem  1  can  now  be  carried  out  exactly 
as  in  Section  2  if  we  make  the  final  observation  that  the 
spectrum  of  U(0)  is  exactly  {x  }.   This  should  be  obvious 


and  follows  from  the  fact  that  on  each  "generalized  eigen- 

j 


space"   E   =©[e   .,U(0)-AI  (for  X  ^  ^^,J)    acts  as 
r      .   r , "]  u 


3 
the  matrix 

x^-X    e     0   .   .   .      0 

0    •    • 

(6.26)  ^  •     ■  Q 

e 

0     .     .     .     0   x  -X 

r 

and  the  inverse  is  obtained  merely  by  inverting  this  matrix. 
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6C.      Extension    to   the   L      Case 

■ oo  

Here  we  would  like  to  make  the  simple  observation  that 
the  difference  equation  theory,  like  the  ODE  theory,  is 
essentially  an  L^  theory  which  produces  AP  solutions  when 
the  coefficients  are  AP. 

In  fact,  if  we  now  agree  that  derivatives  are  to  be 
understood  as  L   instead  of  AP,  the  theory  developed  in 
Sections  3,  4  and  5  is  entirely  valid  in  the  L   case,  with 
only  one  uninportant  exception.   The  single  exception  is 
the  proof  given  in  (3.36)  ff   of  the  convergence  of  the  L^ 

in  the  case 


n 


'   ^n+il-l 


solution  w   to  the  vector  w^  = 
n  n 

I  Yn 

that  the  coefficients  are  not   smooth.   In  fact  if  we 

refer  to  (3.36)  and  (3.37)  we  see  that  (3.37)  need  not 

be  valid  because  the  function  y  which  is  C   (we  take  A  and 

-'  oo 

f  continuous)  need  not  be  uniformly  continuous  on  the 

whole  real  axis .   In  this  case  we  have  not  been  able  to  show 

that  w  ^  w    uniformly  for  all  n  but  only  uniformly  for  nk 
n    n  ■' 

lying  in  compact  intervals. 

To  see  this  let  us  write  w  (k) ,  w  (k)  to  indicate 

n      n 

explicitly  the  dependence  on  k.   Now  if  the  initial  data 
for  w  (k)  is  consistent,  that  is 


(6.27) 


WQ(k) 


k-vO 


r  Yi  1 


I  Yi  J 


for  some  vector  y,  ^  E,  then  it  is  well  known  (see 
Henrici,  p.  244)  that  w  (k)  must  converge   uniformly  on 
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compact  subsets,  to  some  solution  to  the  equation 

(6.28)  y  =  A(t)y  +  f  . 

But  since  every  solution  to  (6.28)  but  one  is  unbounded,  this 
limiting  solution  can  only  be  the  unique  L^  solution  y(t).  Thus 
we  must  show  (6.27)  and  the  proof  of  this  will  in  fact  show 
that  the  vector  y,  in  (6.27)  is  y(0) . 

We  assume  at  least  that  A  is  C   (with  bounded  derivative) 
so  that  the  results  of  Section  5  are  applicable.  We  can  then 
write 

00 

(6.29)  Wq  -  w   =  k   );    G°  g(k,(j-l)k) 

where   g   is    given   in    (3.37),    and  G.    satisfies    the   estimate 

-a   k|n-j| 

(6.30)  llG    II    1  K^    e      -" 

Observe  that  g(k,t)  ->  0  uniformly  for  t  in  compact  intervals.  Let 
T  be  unspecified  for  the  moment  and  rewrite  the  sum  in  (6.29)  as 

(6.31)  w„-w_=  k     I         G°g(k, (j-l)k)+k    I  G°g (k, ( j-1) k) . 

^   ^      I jk|<T  ^  I jki>T   ^ 

The  first  term  ^  0  because  g(k,t)  ->  0  uniformly  on  the  compact 

-a,T 
interval  [-T,T].  The  second  term  is  0(e     )  as  we  can  easily 

see  from  (6.30).  Thus  by  first  choosing  T  sufficiently  large 

to  make  the  second  term  small,  and  then  choosing  k  small  enough 

to  make    the    first   term   small  we   see    that    IIwq-WqII    ->   0     (k   ->   0)    and 

this    establishes    the    convergence   of  w      to  w      uniformly   on 

compact   intervals . 
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6D.   Mean  Value  Properties 


Let     a     e  A 
n 


n 


It   is    shown    in  Corduneanu    (p.    48)    that 


the    limit 
(6.32) 


,    n+N-1 
lim  rr        I         a.    =   m(a^) 


n 


exists    and   that    this    limit   is    uniform   and   independent   of  n, 
which  we   henceforth  set   equal   to    zero.       This    is    of   course 
the   exact   analogy    to   the   mean   value    of    an   AP    function 

T 


(6.33) 


lim 
T->-oo 


f(t)    dt    E   m(f)     . 


0 


Now  if  w   is  our  AP  sequence  solution,  we  would  like  to 
n             ^ 

study  the  behavior  of  m(w  )   as  k  ->  0 .   Define  w^  as 

-^  n  n 


r  y 


(6. 34) 


^n  = 


n+£-l 


n 


where  y  is  the  AP  solution  to  the  ODE 


(6. 35) 


y  =  Ay  +  f 


Clearly 
(6.36) 


m(w^)  = 


f  i"(yn) 


m(y^) 


Now  since  (see  (3.39)) 


(6.37) 


I  w  -w  II 
n   n  oo 


o(l) 
k^l 
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we  must  compare  m(y  )  with  m(y).   We  first  claim  for  any 
AP  function  g  we  will  have 


(6.38) 


"^^^n^  -kW"   ^^^^ 


where  q  =  g(nk).   This  is  obvious  because  (6.38)  holds  if 
^n 

g  is  a  trigonometric  polynomial  and  for  arbitrary  g,  we 
simply  approximate  by  a  sequence  of  trigonometric 
polynomials . 

Now  (6.38),  together  with  (6.37)  and  (6.34),  certainly 
yields 


(6.39) 


m(w  ) 


n   k^O 


f   in(y) 


m(y) 


We   now   concern    ourselves   with    the    rates    of    convergence, 
assuming  smooth    coefficients.       If  A   and    f   are   C  (so    that 

y   is)    we    can    replace    (6.37)    by    (see    (3.41)) 


(6.40) 


w   -w   II       =   O(k^) 
n      n   °° 


where  p  is  the  order  of  accuracy  of  the  scheme.   It 
therefore  follows  that 


(6.411 


m(w  )  = 
n 


n 


m(Y„)  J 


n 


+  0(kP) 


and  it  is  only  necessary  to  study  the  convergence  of  11^(7^^) 


r+1 


to  m(y) .   In  fact  for  any  AP  function  y  e  C    we  have 
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(6.42)  m(y^)    =   m(y)    +    0(k^'^-^) 


Equation    (6.42)    is    an    almost  periodic    analogue    of    a 
theorem   of   Isaacson    and   Keller    [10]    for   periodic    functions 
(p.    340)    and   in    fact   the   proof   is    very    similar,    and   is    a 
simple    consequence   of   Taylor's    theorem.    Letting   t      =   nk 
and   s   s    [0,1]    we   can  write 

(6.43)  y(tj^+sk)    =   y^   +    (sk)y^   +...+   -^^^  y^   +   ©(k"""^^) 

r+1 
where  the  remainder  is  uniform  in  t  since  y     is  AP. 

Integrating  (6.43)  from  s  =  0  to  s  =  1  we  obtain 

^n+1  1 

(6.44)  I   y(t5  dt  =  k  I  y(tj^  +  sk)  ds 

0 

=  k[y^.|y;....  ^-^y;;  +  0(k^^^]. 


^n  0 


If  we    sum    (6.44)     from   n    =    0    to  N ,    divide    by  Nk    and    let  N   ->■   «> 
we   obtain 

(6.45)       m(y)    =    m(y^)    +    |  m(y^)    +    ...         J^Tiyi  '^'^vl')    +    0(k''"^^) 

Now  if  r  =  0,  i.e.  y   is  C    the  terms  involving  ro(y  ) 
through  m(y  )  are  missing  and  we  obtain  (6.42).  For  r  >  0 

we  assume  (6.42)  holds  for  j  =  0,  ...,  r-1.   We  can  then 

1    .    r 
write,  as  y     is  C   , 


(6.46)  0  =  m(y-'-)  =  m(y^)  +  O(k^) 
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and  similar  expressions  for  m(y  )  for  i  =  2,...,r  since 
y  G  c         (we  have  used  the  obvious  fact  that  the  mean 
value  of  a  derivative  is  zero),  to  obtain  (6.42).  If  we  now 
apply  this  to  (6.41)  we  obtain 

m(y) 


(6.47) 


m(w^)  = 


m(y)  ^ 


+  O(k^) 
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6E.   Weakly  Nonlinear  Equations 
Consider  the  equations 
(6.4  8)  y  =  Ay 

(6.49)  y  =  Ay  +  f(t) 

(6.50)  y  =  Ay  +  h(t,y) 

If  any  pair  of  function  spaces  (B,D)  is  admissible  for 
(6.49)  the  contracting  mapping  principle  provides  a  technique 
to  obtain  D-solutions  to  (6.50)  under  certain  conditions  on 
the  nonlinear  term  h(t,y),   the  most  important  being  that 
h  has  a  small  Lipschitz  constant   (see  Hartman  [11], 
Chapter  12)  . 

The  case  of  (A, A)   admissibility  is  particularly  simple. 
We  require  that  h  be  almost  periodic  in  t,  uniformly  for  y 
in  compact  subsets  (see  Hale  [12]  pp.  113  f f ) .   This  simply 
insures  that  h(t,x(t))  will  be  AP  for  any  AP  function  x. 
We  also  require  that 

(6.51)  llh(t,Xj^)  -  h(t,X2)ll^  1  e^Wx^   -   x^W 


fo 


r  any  x,,x  G  E.   The  restriction  on  e^  is  simply 


(6.52)  e  K-,^  <  1 

where  K,  is  the  admissibility  bound  for  (6.49)  (see  (0.3b)) 
The  unique  AP  solution  to  (6.50)  will  simply  be  the  limit 
of  the  iterates 
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(6.53) 


i+1 


is  the  unique  AP  solution  to 


y     =  Ay     +  h(t,y  )  . 

This  is  of  course  an  immediate  consequence  of  the  contracting 
mapping  principle. 

The  same  proof  will  yield  AP  sequence  solutions  to  the 
discretized  version  of  (6.50). 


(6.54) 


w 


n+ 


1  =  ^n^n  +  ^    g(nk,k,w^^^,w^) 


Here    the   nonlinear   term   is 


-1 


(6.55)       g(t,k,w,w)    = 


f    (aj^I-k6^A(t+£k)       [^^h{t,x^_^) 

+      I      6.h(t,x.)] 
j  =  0      ^  J 


f    ^£-l1 


where     w 


x^ 


r    X 


and      w   = 


£-1 


X, 


Observe    that   g   is. 


for    fixed   k,    AP    in    t   uniformly    for  w    and      w   in   compact 
subsets,    and    that   g  will   have    an   0(£,)    Lipschitz   constant 
which  we   call    c,    i.e. 

(6.56)       II  g(t,k,w^,w^)     -    g(t,k,w,,w,) 


'2""2' 


1'    1'    «> 


max  [Ilw2-w^ll   ,  Ilw2-w-|^ll] 


Now  if  K  is  the  admissibility  bound  for  the  linear  inhomogeneous 
version  of  (6.54)  (see  (3.30))  and  \i   =   eK   <    1   then  (6.54) 
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will  have  a  unique  AP  solution  w   ,  which  in  fact  will  be 

n 

the  limit  of  the  iterates  w   .  defined  by 

n;  1  -^ 


(6.57)  w^.Q  =  0 

w   . , ,   is  the  unique  AP  solution  to 

w  ,,  .  ,  =  C  w   .^,  +  kg(nk,k;,w  ,,  .  ,w   .) 
n+l;i+l     n  n;i+l     ^      '  n+l;i   n;i 

This  is  an  immediate  consequence  of  the  contracting  mapping 
principle  and  we  note  for  future  reference  that 

(6.58)  llw  -w    II   =  0(y'")  . 

n   n ;  m  oo 


Now  the  first  thing  to  consider  is  the  convergence  of 

where  y  is  the  solution  to 


r  Yn+£-i 
w   to  w   defined  as 

n     n 

(6.50)  (assuming  of  course  that  (6.52)  holds).   In  fact 

as  was  seen  in  Section  3  w   satisfies  (6.54)  up  to  an  error 

n 

k  o(l)   (0(k^   )    if  everything  is  smooth).  Hence  the 

k->l 

difference      z      =  w   -   w      is    the   unique   AP    solution   to 
n  n        n  ^ 

(^•59)       ^n+1   =   ^n"n   +   ^ [g (nk ,k ,w^^^+    ^n+l'^n^^n^ 

-   g(nk,k,w^^ww    )]    +   k   o(l)     (0(kP'^^)) 
n+i      n 

and  we   have 

(6.60)  11  z    11       <    mII  z    II       +    o(l)     (0(kP)) 

n   0°   —  n   0° 

and   since      y    <    1      this    settles    the   question    of    convergence. 

We   now   examine   the   question   of   the    computability   of 

the   solution    in    the   exponentially    stable    case;    that    is   when 

we    have    constants    K    and   a    such    that 
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(6.61)  IIW^wT-^ll  <  K  e~"^^"  ^^  ,      n  >  j, 

where   W    is  the  fundamental  solution  to  the  homoqeneous 
n 

version  of  (6.54). 

To  do  this  we  use  (6.58)  to  reduce  the  problem  to  a 

study  of  the  effect  of  round-off  errors  and  errors  in  initial 

data  on  the  solution  w     for  a  fixed  m.  This  will  be  a 

n;m 

simple  consequence  of  the  stability  results  of  Section  3 

(see  (3.42)  ff)  for  the  inhomogeneous  equation. 

In  fact  the  results  of  Section  3  show  that  at  the 

first  iteration  we  would  solve  (for  n  >_  0)  for  a  solution 

w   ,  which  is  related  to  w   i  by 
n ;  1  n  ;  1 

(6.62)  w*  ,  -  w   T  +  O(k^)  +  0(e"^'^")  . 

n;  1     n; 1 

We  now  consider  the  effect  of  this  error  on  the  second 

iterate.   Assuming,  for  the  moment,  exact  calculations  and 

* 
exact  initial  data,  we  obtain  a  solution  w  ^,  which  solves 

the  equation 

*  *  *       * 

(6.63)  w.TT^Cw   o+k  g(nk,k,w^,,  i ^w   , )  . 

n+ 1 ; 2     n  n ; 2      ^        n+ 1 ; 1   n ; 1 

* 
If  we  let   z   .  =  w   .  -  w   .  for  i  =  1,2   then  we  obtain 
n ;  1     n ;  1     n  ;  i 

(^•^^)   ^n+l;2  =  Cn^n;2  ^  k [g (nk ,k ,w^^^. ^+  ^^+1; 1 '^n ; l^^n ; 1^ 

-  g(nk,k,w^^^.^,w^.3^)] 

The  last  term  can  be  written  (using  (6.62))  as 
k  O(k^)  +  k  0(e~°''^'^)  and  we  have 


■67- 


n  ,            n     _  1     _  1,  ■ 

(6.65)   z   „  =  k   y   W  wT-^  O(k^)  +  k   y   W  W.  0  (e  ^^) 

2  =  1  ->                                j-1     -■ 

r  -«kn             r        -a  kn 

=  0(k  )  +  0{e  ^       kn)   =   0(k  )  +  0(e  ^      ) 


where  a,  can  be  arbitrarily  close  to  a.  We  have  used 
(6.62),  (6.61)  and  the  trivial  fact  that  t  =  0(e)  for 
any  positive  6 . 

It  is  now  obvious  that  (6.65)  is  preserved  if  we 

* 

include  in  the  equation  for  w   ^  (see  (6.63))  the  effect 

n;  z 

* 

of  roundoff   and  errors  in  the  initial  data  of  w   - 

n;  Z 

(compared  with  the  initial  data  for  w  ,^)  and  it  is  also 
obvious  that  this  process  can  be  repeated  for  m  itera- 
tions, for  any  fixed  finite  m,  and  this  establishes  the 
computability  of  the  solution  w  . 
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Part  II.   Extensions  to  Partial  Differential  Equations 


We  will  now  use  the  techniques  of  Massera  and  Schaffer, 
in  the  simple  exponentially  stable  case,  to  obtain  AP  solu- 
tions to  a  certain  class  of  inhomogeneous  partial  differential 
equations  (PDE's)  to  which  one  can  apply  ODE  formalism  by 
use  of  the  Hi  lie  Yoshida  Theorem.   The  homogeneous  equation 
has  been  extensively  studied  by  Krein  [13],  to  which  we  will 
refer  often.   A  more  succinct  study  of  the  Hille-Yoshida 
Tkeorem  can  be  found  in  the  appendix  of  Lax  and  Phillips  [14]. 

Throughout  Part  II   E  will  denote  an  infinite  dimensional 
Banach  space. 
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7.   Existence  Theorems 
Consider  the  linear  equation 

(7.1)  y  =  By 

where   B   is    an    unbounded  operator.       The    study   of   the   solutions 
to    (7.1)    has   been    answered  by    the   Hille-Yoshida   Theorem   in 
the    case    that   the    resolvent      R;,  (B)       {-    (Al-B)       )    satisfies 
the   inequality 

(7.2)  llR,  (B)ll     <    i      ,  A    >    0    . 

A  —    A 

A  more  general  formulation  is  given  in  Krein  (in  particular 
Section  2,  Chapter  1).   Krein  shows  that  if: 

(7.3)  (a)    B  is  closed 

(b)  B  is  densely  defined 

(c)  liR^(B)ll  < ,    R.P.  A  >  w  , 

^        (Re  A-  CO)  ^ 

then   there   exists    a   semigroup      Z      ,    t   ^  0    ,    such    that: 

(7.4)  (a)         Z      is    strongly   continuous 

(b)  Zq   =    I 

(c)  II  Z^ll     <   Me'^''^ 

(d)  lim   [Z.-I]e/A  converges  iff  e  e  d„ 
A-^0 

and  in  that  case  it  converges  to  Be 

(e)  if  e  e  D_  then  for  t  >  t„  the  unique 

B  U 

solution  to  (7.1)  such  that  y(tQ)  =  e 

is  simply  y(t)  =  Z     e  (uniqueness  is 

t  tQ 
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shown  in  Theorem  2.7,  page  47) 
(f)    Z,  :  D^  ->  D„  and  commutes  with  B  for  t  >  0  . 

t      15        D  — 

For  a  well  posed  problem  we  obviously  require  oj  <_  0  in  (7.4c)  . 
To  insure  exponential  stability  we  are  going  to  proceed  as 
if   0)  =  0   and  introduce  our  own  damping  term  which  will  be 
allowed  to  depend  on  t.   The  reader  should  note  that  if 
0)  <  0   this  term  is  not  necessary.   It  is  well  known  (see 
Krein,  page  4  3)  that  one  can  always  introduce  an  equivalent 
norm  so  that  the  factor  M  in  {7.4c)  can  be  replaced  by  1, 
however  we  will  eventually  approximate  B  by  a  family  of 
bounded  operators   B,   where 

(7.5)  llzj  =  e  "  II  <  M 

and  we  may  not  be  able  to  introduce  a  norm  so  that  M  can 

be  taken  as  1  for  all  h,  and  thus  we  will  leave  M  unspecified. 

Now  to  make  the  homogeneous  equation  exponentially 
stable  we  introduce  a  damping  term  -6 (t)  and  consider 
the  equation 

(7.6)  y  =  [B  -  6 (t) ]y 

For   conditions    on    6    we    first   let   it  be   a   scalar  AP    function 
such    that 

(7. 7)  R.P.  (m(6) )  >  0  . 

rt         ^ 

dr  6 (r)  \      then 
s         ^ 


This  implies  that  if  g(t,s)  =  exp<- 


(7.8)  |g(t,s)  I  <  K  e  "^^  ^^  ,       t  ^  s  , 
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as  can  easily  be  seen. 

It  is  also  possible  to  permit  6  to  be  an  AP  operator 
fianction  satisfying  certain  commutivity  properties.  Speci- 
fically if  Y(t)  is  the  fundamental  solution  to 

(7.9)  Y  =  -  6Y 
then 

(7.10)  (a)    llY(t)Y~-^(s)ll  <  K  e""^^"^^  ,  t>s 

(b)  6(t),  Y(t),  Y~"'"(t)  commute  with  Z^ 

for  all  t  and  all  r  >  0 

(c)  6  (t)  ,  Y(t),  Y~'''(t)  map  D_  into  itself 

a 

and   commute  with   B. 

These  hypotheses  permit  us  to  essentially  treat  6  as  a 
scalar  in  the  following  analysis.   If  we  define 

(7.11)  g(t,s)  =  Y(t)Y"-'-(s) 

then  the  "fundamental  solution"  to  (7.6)  is 

(7.12)  g(t,s)Z^_g  ,  t  >  s, 
and 

(7.13)  llg(t,s)Z^_gll  1  KM  e~^^^"^^  EKe"°'^^"^^ 

after   a   redefinition   of  K. 

Now    consider    the    inhomogeneous    equation 

(7.14)  y   =     [3-6]y   +    f(t) 

where  f  g  A(E)  or  C^(E).   We  want  to  study  bounded  (for  all  t) 
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solutions  to  (7.14).   First  we  see  that  any  such  solution 
is  unique.   In  fact  if  w(t)  were  a  bounded  solution  to 
(7.6)  we  would  have  (for  t  >_  s) 

(7.15)  w(t)  =  Z^_g  g(t,s)  w(s) 

and  w    E    0    follows    immediately   on    letting   s    ->-   -°o . 

It   is    equally   simple    to   see    that   any   C^   solution   must 

be   given  by 

t 

(7.16)  y(t)    =      J    g(t,s)    Z^_^    f(s)    ds    e   c(f)     . 

— oo 

This  is  in  fact  an  immediate  consequence  of  letting  t_  ^  -°° 

in  the  variation  of  constants  formula 

t 
r 

(7.17)  y(t)  =  g(t,tQ)Z^_^  y(tQ)+    ds  g  (t,s)  Z^_gf  (s) 

which  is  proved  in  Krein  (Theorem  6.1,  page  129)  for  the 
autonomous  case,  and  it  is  a  trivial  matter  to  see  that  the 
same  proof  will  work  if  6  is  time  dependent.  Note  that  (7.17) 
need  not  be  a  solution  to  (7.14),  but  any  such  solution 
must  be  given  by  (7.17). 

We  now  consider  the  operator  C  defined  by  (7.16). 
C  is  a  bounded  operator  mapping  L   ->■  L   .   Boundedness  is 
an  immediate  consequence  of  (7.13)  and  in  fact 

(7.18)  llyll       <   K    (=   -)    llfll 

-'       oo     —  Ct  °° 

The  integrand  is  obviously  continuous  if   f  G  C^  and  the 
only  point  in  question  is  its  measurability  for  a  general 
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f  £  L^  (for  measurability  of  functions  in  a  Banach  space 

see  Hille  and  Phillips  [15]).   However  it  is  shown  by  Krein, 

and  will  be  used  crucially  very  shortly,  that  Z   is  the 

Br  ^ 

strong  limit,  pointwise,  of   e    where   B   is  a  bounded 

"Br 
operator   (Theorem  2.9,  page  48).   Since    e  ^  will  be 

norm  continuous  this  establishes  the  integrand  as  the 

pointwise  limit  of  measurable  functions,  hence  measurable. 

We  will  now  show  that   y   as  defined  in  (7.16)  is  a 

solution  to  (7.14).   We  will  make  use  of  the  above  mentioned 

approximation  theorem  of  Krein  (Theorem  2.9).   It  is  shown 

that  there  exists  a  squence  of  bounded  operators  B   such  that 

(7.19)   (a)    II  (B-B  )ell  >  0   (Ve  G  D„) 

*  n      n->^  B 

B  r 

(b)  II  (Z  -z")ell  — — — >  0   Ve  ,  r  >  0,  where  z'^  =  e  "^ 
*        r   r     n-x»  —  r 

(c)  II  Z^ll  <  M  ,   r  >  0 

(d)  B  ,  Z   (r  >^  0)  satisfy  the  same  commutivity 
relations  with  6(t)  (7.10b,c)  as  does  B  and  Z  . 

We  only  have  to  point  out  that  (7.19d)  follows  from  the 

2 
explicit   definition   of   B      as      -A    I    -    A    R,     (B    )    where    A      ->-  °° 

'^  n        n      n  A  n 

n 

(see  Krein,  page  49).   We  point  out  that  the  operators  B 
cannot  in  general  be  considered  as  spatial  discretizations. 
A  family  of  such  approximations  will  be  assumed  in  the 
following  section. 

Now  consider  the  operators  C  (f)  defined  by 
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(7.20) 


'n   =      j    g(t,s)  z"_g  f(s)  ds  E  C^(f) 


It  follows  immediately  that  the  operators  C   are  uniformly 
bounded  operators  mapping  L   into  C   and  in  fact 


(7.21)  lly  II   <  Kllfll 

where  we  take  the  same  bound  as  in  (7.18).   It  is  also 

immediate  that  y    is  the  unique   C    solution  to 

-'n  ^00 

(7.22)  Y^   =    [B^-6]y^    +    f 

and    y      isAPiffis. 
■'n 

We  are  going  to  show  that   C   ^  C  strongly  in  a  certain 
subspace  of  L^.   Specifically  define   S   to  be  the  siobspace 

spanned  by  the  functions  whose  range  has  compact  closure.  It 

is  well  known  (see  Amerio  [16])  that  ACS.   For  f  G  s  we 
will  have 

(7.23)  "^n^^^  ^  ^^^^ 

To  show  (7.2  3)  note  that  by  the  uniform  boundedness  of  C^^ 
it  is  sufficient  to  show  it  for  a  dense  subset  of  S  and  by 
the  definition  of  S  it  is  sufficient  to  take  f  as 

n 

(7.24)  f  =   I   X.(t)e 

j  =  l   ^     ^ 

where   e.  G  E  and  y.  is  the  characteristic  function  of  a 

measurable  set.   By  linearity  we  only  have  to  work  with 

one  term  in  the  sum  in  (7.24).   Thus  if  I  is  some  measurable 
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set    let      f   =    x^e.       Then  w^(t)    =   y(t)    -   y    (t)    is   explicitly 

oo 

(7.25)  w^(t)    =    J    x(t-r)    g(t,t-r)  [Z^-z"]e      dr    . 

0 

Estimating    (7.25)    we   obtain 

oo 

(7.26)  llw   II       <   K    I    e""'^!!  (Z    -z")ell    dr 

n   oo  —        J  r      r 

0 

and  (7.2  3)  is  an  immediate  consequence  of  the  Lebesgue 

dominated  convergence  theorem. 

Observe    that    (7.2  3)    implies    that   C(f)    is    continuous 
if    f   G    S    and   it    is    AP    if    f    is. 

Now  consider  the  almost  periodic  case.   First  we  point 
out  that  if  6  and  f  are  C   (AP  derivatives)  then  y  is 
also   and  in  fact 

(7.27)  y  =  C(-6y  +  f) 
This  is  trivial  if  we  observe  that 

(7.28)  y   =  C  (-6y   +  f) 

■''n    n   ■'n 

and   if  y    is    defined   to  be    the    right-hand   side    of    (7.27)    we 
have 

(7.29)  y-y   =    [C-C    ]  (-6y)    +    C    (-6  (y-y    ))    +    [C-C    ](f)     , 

n  n  nn  n 

and   clearly 

(7.30)  II  (y-y„)ll      — — — >    0 

-^    -^  n      oo      n->"Oo 

Equation    (7.2  7)    is    now   an    immediate   consequence    of   letting 
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n   ->   00    in    the   equation 
(7.31)  •      y^(t^)    -    y^(t2)    =    |      y^(t)    dt 

We   next   show   that   y   satisfies       an  integrated  version 
of    (7.14).       Specifically 


(7.32)       y(t^)    -   y(t2) 

^2  ^2  ^2 


=    B  dt  y(t)    -  dt    6(t)y(t)  +  dt    f(t)     . 

H  h  ^1 

Since  D_,  is  dense  and  B  is  closed  it  is  sufficient  to  show 
(7.32)  when  f  is  of  the  form 


(7.33)  f  =  e^^^  b  ,  ^  ^  °B' 
Now  observe  that  Bf  is  AP.   It  certainly  follows  that 

(7.34)  By  =  C(Bf)  . 

In  fact  this  holds  whenever  Bf  e  C^^  and  is  a  simple  conse- 
quence of  the  closure  of  B,  approximating  the  integral  in 
(7.16)  by  Riemann  sums  over  finite  intervals  and  then 
noting  that  B  applied  to  each  sum  is  the  Riemann  sum 
approximating  C(Bf). 

Now  the  same  property  certainly  holds  for  B   and  C   i.e. 

(7.35)  B^y^  =  ^n^\'^' 

Finally   since  by  its  construction   B    commutes  with  Z 


-77- 


(Krein,  Theorem  2.9)  we  can  write 


(7.36)  B  y  =  C(B  f) 


Since  y   satisfies  (7.22)  we  have 
■'n 


(7.37)      y^(t,)    -   y^(t2) 


=   B^   J    dt   y^(t)    -    j    dt    6(t)y^(t)+    |   dt    f(t) 
^2  ^2  ^2 


=   B    j    dt  y(t)    +    [\-B]     I    dt   y(t)    +    B^  dt(y^-y) 

^2  ^2 


-   I   dt   6(t)yj^(t)    +   I      dt   f  (t) 


f^2 
where    the    fact   that  dt   y(t)    G   D      is    a   consequence    of 

(7.34)    and   the    same    type   of   closure   argument   as    followed 

that  equation. 

If  we   now    let  n   -»-  oo   in    the    last   equation   in    (7.37)    we 

f^2 

see  that  the  term  B„      dt  (y  -y)  >  0   because 

n    J  ■'n   -'         n-H» 

^1 

(7.38)       B     (y    -y)    =     [C    -C](B    f)    =     [C    -C](Bf)     +     [C    -C]((B    -B)f) 
nn  nn  n  n  n 

and   the    right-hand   side    of    (7.38)    clearly   ^    0    in   the    sup   norm 
as   n   ->-  0°.       It   therefore    follows,    on    letting  n   ->  °°   in    (7.37), 
that    (7.32)    does    in    fact   hold   for    functions   of   the    form    (7.33) 
and   thus    for   all    f   s   a. 

Now   if   6    and    f   are   C       (hence   y   is)    we    can   show   very   easily 
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that  (7.32)  implies  (7.14).  In  fact   since  y  is  C  we  can 


let  h  ^  0  in  the  equation           .^,  ^^^ 

t+h                              r  f 

B       f   ds   y(s)                ds    S(s)y(s)  f(s) 

y(t+h)-y(t)                II  J 

(7.39)       ^ =   H H ^   — h 


where   h   can   be   positive   or  negative,  using    the   closure    of  B. 

Similarly    (7.14)    will  hold   if  we    only    require   Bf   is   AP. 

In    fact   since    in   that   case   By   is   AP,    B    can   be   brought   inside 

the    integral    in    (7.32),    which   can    then   be   differentiated 

directly. 

We  point  out  that  in  both  these  cases  both  y  and  By  are 

also  AP.   For  future  use  we  point  out,  using  the  integral 

2 
equation  (7.32)  applied  to  (7.27),  that  if  6  and  f  are  C 

2 
then  y  is  C   and  y  satisfies 

•        • 

(7.40)  y  =  [B-6]y  -  6y  +  f 

i.e.  (7.14)  can  be  differentiated  formally.  (Note  in  parti- 

2 
cular  that  if  Bf  is  AP ,  B  y  is  AP.)   This  argument  can 

P         •    P 
obviously  be  extended;  namely  if  f  and  6  are  C^    then  y  is  C 

and  we  have 

(7.41)  y^  =  [B-6]y^~   +  f^~   -  (6y)'^    +  6y^    ,   r=l,...,p. 


We  will  now  examine  briefly  the  case  that  f  e  S. 
We  still  have  y   ->  y  but  (7.30)  need  not  hold  since  6y  and  f 
need  not  belong  to  S.  However  we  see  very  easily  that  if 
f  e  L   and 
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X  =    C(g) 
(7.42) 

x„   =  c    (g) 
n  n 

then   we  will  have      x    (t)    -^   x(t)    pointwise    and  boundedly.    In 

fact   if  w    (t)    =   x(t)    -    X    (t)  ,    then 
n  n 


(7.43)         llw    (t)  II     <   K 
n  — 


e   "^11  (Z    -Z^)    g(t-r)ll    dr 


and   the   Lebesgue    dominated   convergence   theorem   is    still 
applicable   and  hence    (7.31)    will   still   yield    (7.27)  ,    if  we 
use    the   bounded   convergence    theorem  and    the    fact   that 

(7.44)  y(t)-y^(t)    =     [C-C^]  (-6y+f)  (t)+    C^  (-6  (y-y^)  )  ( t)    ^^:;^>    0 

Thus  y  has  a  bounded  derivative.   Now  if  we  assume  that  f  is 
of  the  form 

(7.45)  f  =  Xjb  ,  b  G  Dg  , 

then  the  derivation  of  (7.32)  from  (7.37)  and  (7.38)  still 
follows,  if  we  use  the  bounded  convergence  theorem  for  inte- 
grals over  finite  intervals.  Since  finite  sums  of  such  f's 
are  dense  in  S,  (7.32)  holds  for  all  f  s  s  and  if  Bf  e  S  or 
6  and  f  are  C  ,  y  will  still  be  the  unique  C   solution  to 
(7,14).   Thus  the  theory  will  hold  for  a  certain  class  of 
bounded,  non  AP ,  inhomogeneous  terms,  but  is  not  as  simple 
as  the  AP  case  and  in  the  future  we  will  restrict  ourselves 
to  this  case. 

Finally  we  point  out  that  (7.32)  can  be  derived  by 
brute  force  dif ferentation  of  the  formula  (7.16)  but  use  of 
the  approximating  operators  B   permits  a  cleaner  and  more 
straightforward  development. 
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8.  Discretization   in  Space 

The  operators   B   which  proved  so  useful  in  the 
preceding  section  need  not  correspond  to  a  discretization 
in  space.   We  therefore  postulate  the  existence  of  a 
family  of  bounded  operators   B,  ,  where  h  will  generally 
play  the  role  of  a  spatial  grid  size,  such  that: 

(8.1)   (a)    B,  are  defined  and  bounded  for  h  6  (0,hp.] 

(b)  llB,  II  =  O(-i-)  ,  m  >  0 

n       j^ra 

(c)  B,       generate      bounded   semi-groups     (i.e. 

B    t  , 

He   ^      =    zjll     <    M    for    t    >    0) 

(d)  e   e   D^    implies    II  (B-B,  )ell       .     -,    >    0 

B        '^  h  h^-0 

(e)  If   6    is   not    a   scalar    function,    then   the 
commutivity   relations     (7.10)b,c      hold  with 

1- 

B  and  Z   replaced  by  B,  and  Z 
r   '^       -^   h      r 

(f)  If  e  s  D^  (or  any  dense  manifold  in  D^ 

invariant  under  Z   for  r  >  0)  we  have 

r        — 

IIb,  Z  ell  _<  K(I,e)  for  r  e  I,  where  I  is  any 
compact  interval  on  the   nonnegative  real  axis. 

We  point  out  that  (8. If)  is  designed  to  insure  that 

Z   ->•  Z   strongly,  for  all  r  >_  0.   This  is  certainly  reasonable 

since  B,  Z  e  converges  to  BZ  e.   It  is  of  course  trivial  if 
h  r        ^        r 

B,  commutes  with  Z  .   Also  note  that  the  convergence  in  (8.  Id) 
h  r 

does  not  deal  with  the  order  of  the  approximation  of  B,  to  B. 
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This  will  be  dealt  with  by  a  subsequent  assumption. 

Before  proceeding  further,  we  indicate  a  simple  example 
of  the  case  we  are  considering. 


(8.2)   (a)    E  =  L2  (0,27t) 

(b)    B  =  T—   (with  periodic  boundary  conditions) 

0  X 

,  .    _,  _  u(x+h)  -u(x-h) 
(^)    V  2h 


In  this  example  (8.  If)  is  unnecessary  as  B,  in  fact  commutes 

with  Z  . 
r 

Now  associated  with  B,  we  have  the  equations 

(8.3)  y  =  [Bj^-6]y 

(8.4)  y  =  [Bj^-6]y  +  f 

The  homogeneous  equation  (8.4)  has  as  its  fundamental  solution 

B^t 
y,  (t)  =  g(t,s)  e    ,  and  is  exponentially  stable,  i.e.  there 

exist      constants      K    and   a      independent   of      h      such    that 

(8.5)  IIY,   (t)Y"-'-(s)ll     <   K   e"°'^^"^^     ,  t    >    s. 

n  h  —  ~ 


Associated  with  (8.4)  we  form  the  operator  C,  where 

t 
(8.6)      y^  =  C^U)    =      J  g(t,s)  zj_g  f(s)  ds 


for   any  AP    f,    y,     is    the    unique   AP    solution    to    (8.4). 
Equation    (8.5)    shows    that    the   operators      C,     are   bounded 
uniformly   in   h,    i.e. 

(8.7)  lly^ll       <    Kllfll 

where   K    is    independent   of  h. 
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Now  our  first  task  is  to  show  that  C,     ^  C   strongly. 
Using  the  results  of  Section  7  (see  (7,23)  ff)   it  is  suffi- 
cient  to  show  that   Z   ->  Z   strongly,  for  r  >_  0.   This  will 
be  an  easy  consequence  of  (8. If).   In  fact  if  e  ^    D^    (or 
the  manifold  described  in  (8. If))  and  if  y  =  Z  e  , 
y,  =  Z  e  ,  then  if  w  =  y~yv^  we  will  have 

=  Bj^Wj^  +  [B-B^]y 

(8.8)  , 

w  (0)   =   0  . 


From  (8.8)  we  can  write 

t 

(8.9)  w^(t)  =  I  ds  Z^_g [B-B^lZ^e 

0 

and   the    fact   that   II  w    (t)ll    -^    0 ,    in    fact   uniformly   on    compact 
intervals,    is    an   immediate    consequence    of    (8. If)    and   the 
Lebesgue   bounded   convergence    theorem.      The    result    for    any 
e   e   E      follows    from   the    uniform  boundedness    of    the 
operators      Z        and    the   denseness    of    the   manifold   of    (8. If). 
Now   the    fact   that 

(8.10)  II  (y-yj^)IL  ^^7o->   0 

for  all   f  G  A   (in  fact  for  all   f  e  S)  is  not  any  improvement 
of  the  corresponding  result  for  the  operators  B   of  Section  7. 
We  want  to  show  that  the  difference  in  (8.10)  approaches  zero 
as  fast  as  some  power  of  h,  and  to  do  this  we  must  require 
that  B,   approximate  B  up  to  a  certain  degree  of  accuracy. 
Specifically  we  suppose  that  there  exists  an  operator  L 
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(unbounded)    with      D^    c    d„    and   such   that   if   e   s    d,    we   will 

Li     —         O  Li 

have 

(8.11)  II  (B-B,  )ell    <    C  h^llLell       where      j    >    0 

h    — 

and  C  and  j  are  independent  of  h  and  e. 

If  we  have  Lf(t)  AP  we  would  like  to  show  that  (8.11) 
implies 

(8.12)  II  (y-yj^)ll^  =  O(h^)  . 

Unfortunately  (8.12)  requires  a  further  assumption;  namely 


(8.13) 


(a)  Z^:  D^  -  D^ 


(b)    If  Lf  is  AP  then  llLZ^f(t)ll   will  be 


uniformly  bounded  for  r  >_  0  and  all  t. 

Assumption ( 8. 13)  is  of  course  trivial  if  L  commutes  with  Z  . 

It  is  now  a  simple  matter  to  show  (8.12).   In  fact 
w,  =  y~yu  is  the  unique  AP  solution  to 

(8.14)  w^  =  [Bj^-6]w^  +  (B-B^)y 
whence  by  (8.7)  we  have 

(8.15)  liw,  I!   =  0(11  [B-B,  ]yll  ) 

and  so  it  is  only  necessary  to  study  (B-B,  )y.   But  now  if 
Bf  G  A,  which  will  generally  follow  if  Lf  G  A  since  L  will 
be  a  higher  order  operator  than  B,  we  can  write 
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(8.16)  [B-Bj^]y   =    J    dr    g  ( t,  t-r)  [B-Bj^]  Z^f  (t-r) 

0 

and   clearly 


(8.17)  II  [B-B,   ]yll       =   O(h^) 


is  an  immediate  consequence  of  (8.13) 
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9 .  Discretization  in  Time  by  the  Multistep  Linear  Method 

We  will  now  apply  the  theory  developed  in   Section  3 
to  the  spatially  discretized  equation 


(9.1) 


Yh   =    tV^^  ^  f 


Since  B,  is  bounded  the  theory  is  immediately  applicable  and 
we  can  assert  that  for  k  sufficiently  small  (but  depending 
on  h)  there  will  be  an  AP  sequence  w  (k,h)  which  will  converge 
to  w  (k,h)  defined  by 


,-  yj^((n+£-l)k) 


(9.2) 


w^(k,h)  = 


^   yj^(nk) 
and   if  everything   is    smooth    in    t  we  will   have 


(9.3) 


II  w   -w    II       =   0,   (k^) 
n      n   °°  h 


where  the  bound  in  (9.3)  may  depend  on  h.   It  is  only  left 
to  show  that  the  bound  in  (9.3)  is  independent  of  h  and  to 
discuss  the  relationship  of  k  to  h.  We  will   show  that  what 
is  required  is 


(9.4) 


k  + 


2m 


0(1) 


i.e.  there  exists  a  constant   c   independent  of  k  and  h, 
such  that  if  k  +  k/h    <  c   the  difference  equation  corres- 
ponding to  (9.1)  will  be  exponentially  stable  with  constants 
independent  of  k  and  h. 
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The  quadratic  dependence  on  h  in  (9.4)  is  restrictive 
and  will  be  removed  in  the  following  section. 

Throughout  this  section  we  will  use  the  notation  of 
Sections  2  and  3.   We  will  also  take  6  to  be  C  . 

We  now  fix  h  and  consider  the  homogeneous  difference 
equation  on  E„  (see  (3.7)), 

(9.5)  w  ^,  =  U(kB,-  k6  )  +  k  0(l)w   . 

n+1       h    n  n 

It  is  clear  that  this  will  hold  uniformly  in  k  and  h  if 
k  +  k/h   is  sufficiently  small,  i.e.  if  we  are  in  a  region 

(9.6)  k  +  ^  =  0(1)  . 

h 

Of  course  restrictions  of  the  form  (9.6)  are  included  in  (9.4) 
Also  note  from  the  explicit  definition  of  the  full  linear 
term  in  (9.5)  (see  (2.11))  the  last  term  in  (9.5)  will  be 
uniform  in  k  and  h.   For  the  rest  of  this  section  this  will 
be  understood  for  all  "O"  signs  unless  stated  otherwise. 

If  W  (k,h)  is  the  fundamental  solution  to  (9.5)  it  is 

n 

clear  that  all  we  must  show  is  that  in  a  region  of  the  form 
(9.4)  we  will  have 

_i         -a  k(n-j) 

(9.  7)  IIW  W.  II  <  K^  e   -^        ,      n  >  j  , 

n  j    —  1  — 

where  the  constants  K^  and  a,  are  independent  of  k  and  h. 
Following  the  procedure  of  Section  3  we  let 


(9.8)  v  =  T(kB,  -  k6,  )w^ 

n       n     n  n 
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and  we  then  have  in  a  region  of  the  form  (9.6) 


(9.9) 


v„,,  =  L(kB   -  k6  )v^  +  k  0(1)  v^ 
n+l       n     n   n  n 


Using   the    structure   of   the  block   diagonal   operator      L(z) 
derived   in  Theorem   1     we   can  write    this    as 


(9.10)       V 


n+l 


(    [I+k(B    -6    )  0 

n      n 

0  B 


n 


V      +   kO  k   +   -T—  V 
n  ,  2m      n 

h 


Here  we   have    used   the   notation   B      for   the    lower  block    (see 

n 

(2.26))  so  as  not  to  confuse  it  with  the  operator  B,  . 
Observe  that  in  (9.10)  we  have  the  first  appearance  of 
the  term  k/h   .   This  is  because  in  order  to  apply  Theorem  2 
we  must  separate  out  a  factor  of  k  in  the  perturbing  term. 
How  it  certainly  follows  that  in  a  region  of  the  type 

(9.6)  we  have 


(9.11) 


n-1 

!   '.   1   —   1 


n-D 


n  >  j 


where   K,  and  6  are  independent  of  h  and  k.   We  therefore 
only  have  to  consider  the  upper  block  in  the  leading  term 
of  (9.10)  and  in  order  to  apply  Theorem  2  we  are  going  to 
compare  this  with   Y   ,  =  Y,  (nk)   where   Y,  is  the  fundamental 
solution  to 


(9.12) 


^h  =  ^\-'^\ 


Observe    that   the   estimate 

-1 


(9.13) 


Y      ,  Y    T   ,  II     <    K   e 
n ;  h      :) ;  h 


-ak (n- j) 


n    >    j 


holds,  where  K  and  a  are  independent  of  h. 

We  will  show  that   Y   ,  satisfies  the  difference 

n;h 

equation 


(5-1^)    ^n-.l;h  =  fl^^^n^'^'^^l^^-h 

=  [I+k(B^-6^)]Y^^.^  +  kO(k  +  ^)Y^^^    . 

Observe  that  the  last  term  in  (9.14)  is  of  the  same  order 
as  the  perturbation  in  (9.10). 

Equation  (9.14)  follows  easily  from  the  same  argument 
as  in  Section  3  (see  (3.13)  and  (3.14))  if  we  simply  observe 
that 

(9.15)  (B^  -6)  =  -  6 

and  that  (3.15)  holds  in  this  case  for  t  >^  s  by  (8.5). 
We  can  now  apply  Theorem  2   comparing  the  equation 

(9-1^)  ^n+l;h  =  ^^^^^\-'n^^\;h 

with    (9.14).      We   simply   note    that   the   proof  of   Theorem   2 
indicates    that    the    terms      £„    ,    a,    and  K,     (see    the    statement 
of   Theorem   2)    depend  only   on    the   constants    K   and   a   in    (9.13) 
and    are    thus    independent   of  both   k    and   h.    Another   applica- 
tion  of   Theorem   2   comparing   the   equation 


(9-1^)  Vl 


f    l+k{B^-6J       0 

0  B 

n 


V 
n 


with    (9.10),    together  with    the   boundedness    of    the    operators 
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T(k{B,-6  )),  T~''"(k(B,-6  ))  in  the  region  described  by  (9.6) 
n   n  n   n  ■' 

will  yield  (9.7)  with  the  restriction  k  +  k/h    sufficiently 
small  (i.e.  (9.4) ) . 

If  we  now  refer  to  the  proof  of  admissibility  in  the 
stable  case  given  by  (3.31)  ff  we  see  immediately  that  (9.7) 
implies   (L^,L^)  admissibility  for  the  inhomogeneous  version 
(9.5),  with  an  admissibility  bound  independent  of  k  and  h. 
It  is  only  necessary  to  show  that  the  convergence  of  w   to  w 
is  independent  of  h  (see  (9.3)).   This  will  follow  from  (3.40) 
if  we  can  show  Hyp   11^  will  be  bounded  uniformly  in  h, 
assuming  that  6  and  f  are  smooth. 

In  fact  if  6  and  f  have  p+1  derivatives  the  differenti- 
ability properties  of  y,  are  the  same  as  those  of  y   introduced 
in  Section  7.   In  addition,  the  formula   given  in  (7.28)  can 
obviously  be  extended  (for  the  operators   C,  )  to  yield 

(9.18)    yf  1  =  C^(fP^l  +  (-6y^)P-*-^+  6yP^S 

(note  that  the  argument  of  C,  in  (9.18)  does  not  involve  y^   ) 

and  the  uniform  boundedness  of  the  operators  C,  immediately 

yields  the  uniform  boundedness  of  the  derivatives  of  y,  . 

f   y((n+£-l)k  " 


If  we  now  define  w   = 

n 


,  then  (9.3) 


y(nk) 
(uniformly  in  h)  together  with  (8.12)  implies 


C9.19)      ""n"%"oo  =   0(kP+h^)  =  0(h^"^P+h^) 
where  we  have  used  (9.4)  for  small  k. 


-90- 


10.   The  Lax-Wendroff  Scheme 

The  term  k/h    in  (9.4)  can  be  improved  for  operators 
such  that  the  Lax-Wendroff  scheme  is  stable  for  the  homogeneous 
equation 

(10. li  w  =  Bw 

We  will  then  be  able  to  obtain  AP  sequence  solutions  under 
the  more  favorable  restriction 

(10.2)  k  +  ^  =  0(1)  . 

h 

3 
First  of  all  we  require  that  6  and  f  are  C  .   As  shown 

3 
in  Section  7  this  implies  that  y  is  C  .   It  implies  further 

that  y  satisfies 

(10.3)  y   =    [B-5]y    -    6y   +    f 
i.e.    that    the  basic   equation 

(10.4)  y    =     [B-6]y    +    f 

can   be    differentiated    formally. 

We    can   now   apply    the   Lax-Wendroff   scheme    to    (10.4). 
Proceeding    formally  we  write,    assuming   Bf   is   AP, 

2 
2 


(10.5)       y(t+k)    =   y(t)    +    ky(t)    +    5-  y  (t) 


k^  2 

=   y(t)    +   k[B-6]y(t)    +   J-    [B-6]    y{t) 

k^         • 
+   ^    [-6y   +(B-6)f   +f]    +   kf(t)    = 
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k^   2 
=  [I+kB  +  y-  B^]y(t)  -  k6y 

+  k[-5-  6^-  y  6  -  k6B]y+  k[f+  y(B-5)f+  |  f] 


We  have  written  this  as  an  inhomogeneous  difference  equation 
in  order  to  apply  techniques  already  developed. 

Now  in  the  expression  for  y  we  approximate  B  by  B, 

2      2 
and  in  the  expression  for  y   approximate  B   by  B  ,  ,  for  a 

certain  constant   c  >  0   to  be  determined  (c  =  1/2  for  the 

example  described  in  (8.2)).   Of  course  we  will  have  to 

2  2 

assume  that   B,  approximates  B   in  the  same  sense  that  B, 

approximates  B,  but  we  will  leave  the  precise  hypothesis 

for  later. 

Applying  these  approximations  to  (10.5)  and  letting  w  (k,h) 

stand  for  the  dependent  variable  we  obtain  the  difference 

equation 

(10.6)  w  ^T  (k,h)  =  U(k,h)w  -  k6  w  +  kR(n,k,h)w^+  ki 

n+1  n     n  n  n     n 

where 

k^   2 

(10.7)  (a)    U(k,h)  =  I  +  kB^  +  2-  ^ch 

(^)    V  ^n  -^  I  t^ch  -  ^n^^n  "^1  ^n 

(c)    R(n,k,h)  =  I  6^  -  I  6^  -  k6^B^^  =  O  (k+  ^)  . 

Now  U(k,h)  represents  an  approximation  to  (10.1).  We 
choose  c  so  that  for 

(10.8)  -~   =  0(1) 

h 
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this    discretization   is    stable,    i.e. 

(10.9)  IIU(k,h)^ll     <    H    ,  r   >_  0. 

We    are    going   to   use   Theorem   2    to   show   that   if   liR(n,k,h)ll 
(i.e.    k   +   k/h    )    is    sufficiently    small,    then    the   homogeneous 
version   of    (10.6)    is   exponentially   stable;    that    is    if 
W    (k,h)    is    the    fundamental   solution    to   the   homogeneous 
equation,    then    there   exists    constants      K,    and   a,    independent 
of  k    and  h   such    that 

_-,  -a,k(n-j) 

(10.10)  IIW  W.    II     <    K,    e  ,  n    >    j, 

n    3        —     1  '  —  -^ 

where  for  simplicity  we  have  suppressed  the  dependence  of  W 
on  k  and  h.   Once  we  have  obtained  (10.10),   (L  ,L  )  (and 
(A^fA^))  admissibility  will  follow  in  the  usual  manner. 

To  obtain  (10.10)  we  observe  that  since  we  can  allow 
k  +  k/h   to  be  as  small  as  required.  Theorem  2  permits  us 
to  neglect  the  perturbing  term  R.   Also  by  neglecting  a 
perturbation  of  the  same  order  as  kR  we  can  work  with  the 
simpler  equation 

(10.11)  w^,,  =  [I-k6^]  U(kh)w^  . 

n+1        n        n 

Now  consider  the  function   g(t,s)  defined  in  (7.11). 
Using  the  same  argument  given  in  Section  3  (see  (3.16)),  we 
see  that   g   =  g(nk,0)  will  satisfy  the  equation 

(10.12)  g^+i  =  f^+^^n^^n  =  [I-k6^+k2o(l)]g^ 

where  the  "O "  sign  in  (10.12)  is  independent  of  n.   Thus  by 
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again  neglecting  a  perturbation  of  the  form  kO  (k+k/h"^) 
we  obtain  finally  the  equation 


(10.13)  w  ^,  =  [I+kD  ]  U(k,h)w   . 

n+1        n      '    n 

But  the  fundamental  solution  to  (10.13)  is  simply 

_-l 

(10.14)  W^W.    =  g(nk,  jk)U(k,h)^~^ 

as    one   can   verify   immediately,    using    the   comrautivity 

properties    of      g     with      B,        (see    (8.1e))    in   the   case   that 

n  h 

6    is    an   operator.      Now   the   estimate 

_-l  -ak  (n- j) 

(10.15)  IIW  W.    II     <    K   e  ',  n    >    j, 

n    j        —  —   -" 

follows  from  (7.10a)  and  (10.9)  and  we  can  now  apply 
Theorem  2,    noting,  again  that  the  terms   K,  ,  a,,  e^  depend 
only  on  the  terms  K  and  a  (see  the  hypothesis  of  Theorem  2) , 
to  obtain  (10. 10) . 

oo    oo  ^     ^ 

We  now  have  (L  ,L  )  (or  (A  ,A  ) )  admissibility  for 
n   n        n   n 

(10.6)  and  it  is  only  necessary  to  study  the  convergence  in 

the  case  that  the  inhomogeneous  term  is  given  by  (10.7b). 

2  2 

To  do  this  we  must  now  assume  that  B,  approximates  B  ;  namely 

that  there  exists  an  operator  L  such  that  if  e  s  d   n  d  ^ 

L  ^ 

then 

(10.16)  II  (B^-B^)ell     <   C^h^llLell 

n  —2. 

We  also  assume  that  L  satisfies  the  same  assumption  as  L 
(see  (8.13))  so  that  if  Lf  and  Lf  are  AP  we  can  conclude 
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(using   the    same    argument   as    in    (8.16)    and    (8.17)) 


(10.17)        (a)         II  (B-Bj^)yll^   =   O(h^) 


(b)  II  (B^-B^)yll^   =    O(h^) 


If     w    (k,h)    is    the    unique   AP   solution    to    (10.6),    we    can   now 
n  ^ 

show    that     w      converges    to   y      uniformly   in   n. 

3 
In    fact   since   y   is   C    ,    y      satisfies     (10.5)    up   to    an   error 

3 
0(k    ).    If  we    then   make   use   of    (10.17)    we   see   that   y^   will 

satisfy    (10.6)    with   an   error   kO(h-'+k    ),    where    this    is    uniform 

in   n.    Then    as    the    admissibility   bound   of    (10.6)    is    independent 

of   h     (by    (10.10)),    we    can    conclude 


(10.18)  llw  -y   II      =  0(h^+k^)    =  0(h^+h^"^) 


n      n   0° 

where  we  have  used  (10.2). 

For  the  example  described  in  (8.2)  we  will  have 

3  4 

9         ~     3  1 

(10.19)  L  =  J   and  L  =  — ^  ,c  =  2'"^=l'j  =  2, 

8x  8x 

(periodic  boundary  conditions) 

Finally  we  would  like  to  point  out  that  in  the  term  f^ 
in  (10.7b)  one  can  replace   B^^f^  by  Bf,  without  changing 
the  order  of  the  error.   This  is  simply  because  Lf  is  AP , 
hence 

(10.20)  li  (B-B^)fll^  =  O(h^)  . 
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11.   Linear  Perturbations 
We  now  consider  the  perturbed  system 

(11.1)  y  =  [B-6]y  +  D(t)y  +  f(t) 

Here  D(t)  is  a  bounded  operator  for  each  t  and  is  almost 
periodic  as  an  operator  function.   We  also  require 

(11.2)  llDll   <   e. 

oo   (J 

where  Cq  is  to  be  specified. 

The  first  step  in  finding  AP  solutions  to  (11.1)  is 
to  solve  the  equation 

(11.3)  y  =  C(Dy)  +  C(f)  . 

This  equation  can  be  solved  uniquely  by  the  contracting  map- 
ping principle  provided  y  =  e-K  <  1,  where   K  is  the  norm 
of  C  (see  (7.18)).   The  contracting  mapping  principle  also 
shows  that  the  assignment  of  y  from  f  is  a  bounded  operator 
and  in  fact 

(11.4)  llyll        <   ^   II  fll       =   KJIfll 

^        '  Jf   00  —  1-y       oo       1      00 

where  K,  depends  only  on  e„  and  is  independent  of  D. 

We  want  to  show  that  y  satisfies  (11.1),  and  if  we 
assume  that  D,  6  and  f  are  C  ,  the  results  of  Section  7 
show  that  it  is  sufficient  to  show  that  y  is  C  .  This  will 
follow  if  we  observe  that  y  will  be  the  limit  of  iterates 
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y.    defined    as 

(a)  y      E    0 

(11.5)  " 

(b)  y^^^   =    C(Dy^)    +    C(f)     . 

We  will   then  have    for  each   i 

(11.6)  y^^^   =   C(Dy^)    +   C(-6y^    +   Dy^    +    f) 

If  we  then  define  y  as  the  unique  solution  to 

(11.7)  y  =  C(Dy)  +  C(-6y  +  Dy  +  f) 
we  will  have 

(11.8)  y^  -^  y. 

To  see  (11.8)  note  that  H y • H ^  is  bounded  uniformly 

in  i  (from  (11.6)  and  y  <  1)  and  if  g.  is  defined  as 

1 

(11.9)  g^   =    sup   lly-y^ll 

r>i 


Th 


en      {g.}    form  a      nonincreasing   sequence    satisfying 


(11.10)  g^^^  1  u  g^  +  o(i) 

i-voo 

and  (11.10)  implies  g.  -^  0  which,  in  turn,  implies  (11.8), 
The  fact  that  y  =  y   follows  exactly  as  in  Section  7  (see 
(7.27)  ff)  .   Observe  that  since  y  is  AP ,  By  is  also  AP 
(from  (11.1) ) . 

Now  consider  the  operators  C,  defined  by  (8.6).  The 
equation 

(11.11)  y^  =  Cj^(Dyj^)  +  C^(f) 
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can   be    solved  exactly   as     (11.3)    can.       (For   simplicity  we 
suppose    that   HC,  II    =    Hcll    so   that   we   have   the   same    restric- 
tion  on    £«• ) 

Our  first  task  is  to  show  that  y,  -^  y.  This  will  be  a 
simple  consequence  of  the  strong  convergence  of  C,  to  C. 
In  fact  if  w,  =  y  -  y,  we  have 

(11.12)  w^   =    [C-Cj^]  (Dy)    +    C^(DWj^)    +     [C-C^]  (f) 
whence 

(11.13)  llw,  II       <    ijllw.  II       +    o(l) 

h   °°  —  h   °°        ,     -, 

h^O 

and   this    certainly   implies      llw,  II      >    0. 

"        h->0 

Now   of   course    this   would  not  be  useful   unless   we    could 
prove 

(11.14)  llw^ll       =   O(h^) 

h  °° 

This  estimate  can  be  shown  provided  we  make  certain 
assumptions  on  the  interaction  of  D  and  the  operator  L 
defined  in  (8.11).   We  first  observe  that  w,  satisfies 
the  equation 

(11.15)  Wj^  =  Cj^(DWj^)  +  Cj^([B-B^]y) 

and  using  the  analogue  of  (11.4)  for  the  equation  (11.11), 
together  with  the  uniform  boundedness  of  the  operators  C,  / 
we  see 

(11.16)  llw^ll^  =  0(11  [B-Bj^]yll^)  , 
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and  it  is  sufficient  to  show 

(11.17)  II  [B-B^]yll^  =  O(h^)  . 

To  obtain  (11.17)  we  make  an  assumption  on  the  operator 
L.  We  will  call  it  Assumption  I  as  we  will  later  have  to  make 
the  same  type  of  assumption  for  other  operators. 

Assumption  I:   Suppose  there  exists  a  sequence  of 
operators   L»  =  I,  L, / . . . ,L  =  L  such  that 

(11.18)  (a)         L.    are   closed 

(b)  E   =   D^       3   D^       3    ...    3   D^ 

^0   ~     ^1  ~  "        n 

(c)  L.    commutes   with    Z.     for   t    >    0    (i.e.    Z.  :    D^    ->D^ 

*    '  1  t  —  t        L.       L. 

and        L.Z.e   =   Z.L.e    for  e  e  d^.    ) 
It  t   1  i-i  ■ 

(d)  If  e  G   D^      then   D(t)e  e   d^      and 

?       a 

L   D(t)e    =       I      E  .  (t)L.e 

^  j=0       ^  ^ 

£ 
where    the    operators    E.(t)    are   bounded   and  AP, 

and   eJ   =    D. 

Assumption  I  is  clearly  motivated  by  the  example  where 
D(t)  is  multiplication  by  some  function  d(t,x)  (x  is  the 
spatial  variable) . 

Now  Assumption  I  implies  that  if  L^^f  is  AP  then  L^^y  is  AP , 
for  £  =  l,...,n.   First  note  that  if  z  =  C(f)  then 

(11.19)  h^z   =  C(L^f)  . 

This  is  a  consequence  of  (11. 18a, c)  using  the  same  argument 

as  was  used  for  B  (see  (7.34)  f  f)  .   Now  if  we  define  Z£.j_=  L^^y^ 

where  the  y.  are  defined  in  (11.5);  then  it  certainly  follows 
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from  (ll.lSd)  that   z   .   is  AP  for  i  and  for  £  =  0,...,n. 

A/  f  -L 

We  will  now   show  that  the  sequence  {z„  .}  converges  as  i-^ 
To  see  this  first  set  £  =  1 .   We  then  have 

(11.20)  ^j^.^+i   =    C(Dz^.^)  +  C(Ej(t)y^)  +  C(f) 
and  z,  .  ^  z,  ,  where  z   is  the  solution  to 

(11.21)  z^  =  C(Dz^)  +  C(Ejy)  +  C(f)  . 

In  fact  w^  .  =  z^  -  z^  .  satisfies  the  equation 

(11.21)      ^i,i+i  "  C(Dw^^^)  +  C(Ej(y-y^)) 
and 

(11.22)        "w^,iiL  =  "^r^i,i"a>  -i^^>  0 

follows  from  the  same  argument  used  in  proving  y  is  C 
(see  (11.9)  f f )  .   Finally  by  the  closure  of  L-j^  we  must 
have  z,  -   L.y.   It  is  clear  that  a  simple  induction,  using 
the  same  proof,  will  establish  that 

(11.23)  L^y^   i^^  >  L^y  ,  i   =    l,...,n. 

In  particular  L  (=  L  )y   is  AP  and  (11.17),  hence  (11.14) 
is  valid. 

Finally  we  point  out  that  if  we  assume  a  chain  of 
operators,  as  in  Assumption  I,  for  the  operator  B,  we  will 
then  have  BDy   AP.    We  will  also  have  to  make  such  an 
assumption  for  the  operator  L  (see  (10.16)),  when  applying 
the  Lax-Wendroff  scheme  to  (11.1),  as  will  be  discussed 
shortly. 
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We  will  now  apply  to  (11.1)  the  two  methods  of  time 
discretization   that  have  already  been  described.   We  are 
interested  in  the  existence  of  AP  sequence  solutions,  or 
equivalently  exponential  stability  for  the  homogeneous 
difference  equation,  with  k  and  h  restricted  according  to 
relations  (9.4)  and  (10.2).   As  one  might  expect   this  will 
be  a  simple  consequence  of  Theorem  2,  provided  £„  is 
sufficiently  small. 

We  consider  first  the  linear   multistep  scheme  as 

discussed  in  Section  9.   Since  the  linear  term  C   of  the 

n 

difference  scheme  (see  (2.11))  is  a  smooth  function  of  its 

£+1  arguments,  it  follows  that  if  k  +  k/h   is  sufficiently 

small  (i.e.  IlkB,  II  is  small)  then  the  difference  equation 
h 

for  (11.1)  can  be  written 

(11.24)      w  ^T  =  C*w   +  k  0{e^}w^    +  ^^  f„ 
n+1     n  n         On       n 


where   C    is  the  linear  term  corresponding  to  the  unperturbed 
version  of  (11.1),  i.e. 


(11.25)  y  =  [B^-6]y  +  f 

Since  we  have  shown  in  Section  9  that  in  a  region  of  the  form 

(11.26)  k  +  -|j^  -  0(1) 

h 

the  unperturbed  homogeneous  difference  equation 

(11-27)  w^^^  =  C*  w^ 

is  exponentially  stable,  with  constants  K^^  and  a-,^  independent 
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of  k  and  h  (see  (9.7)  ff)  it  will  follow  from  Theorem  2 
(see  the  remark  following  (9.16))  that  in  the  same  region 
the  homogeneous  version  of  (11.24)  will  be  exponentially 
stable  if  e^   is  small  enough.   It  then  follows  easily 
that  the  inhomogeneous  equation  (11.24)  will  have  (L^,L^) 
(or  (A  ,A  ))  admissible,  with  an  admissibility  bound 
independent  of  k  and  h. 

We    can   thus    obtain  AP    sequence   solutions    to    (11,24) 
and,    as    shown   in   Section    9    (see    (9.18)    f f ) ,    we  will   obtain 
convergence    of   O(k^)    provided        II  y^      11^      is   bounded   in   h. 
However,    as    one   can   trivially   see,    if  D,    6,    and   f    are   C^ 
then    the    same   argument   used   to   show  y   is   C       (see    (11.6)    ff) 
can  be    repeated  p+1   times    and   applied   to   the   operators   C, 
as   well   as   C    to   show   that  y   and  y,     are   C^        and   in    fact    (11.7) 
can  be    generalized   to    (compare  with    (9.18)) 

(11.28)      y^   =   C(Dy^)    +   C([(D-6)y]^    -    (D   -    6)y^   +    f^)     , 

i    =    1, . . . ,p+l, 

with  a  similar  formula  for  y,  and  C,  .   This  certainly  shows 

h      n 

that   lly,  II   will  be  bounded  in  h  for  i  =  l,...,p+l   and 
thus  that  we  can  obtain  a  complete  extension  of  the  theory 
of  Section  9.   We  note  in  passing,  that  the  same  proof 
which  led  to  (11.13)  will  show  that   Yh  "^  ^  * 

The  extension  of  the  Lax-Wendroff  scheme  will  be  just 
as  simple.   In  fact  (under  the  assumption  that  D,  6  and  f 
are  C  )  the  difference  equation  for  (11.1)  will  have  the 
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same    linear   terra   as    the   equation    for    (11.25)  ,     (10.6) 
except   for   a  perturbation 

(11.29)  k[D^    +    I    ((B^^    -    6^)D^    +    D^(B^^    -    6^)    +    D^    +    D^) ] 

as  one  can  easily  see  by  carrying  out  the  expansion  of 
Section  10  for  the  equation  (11.1).   It  follows  that  if 
k/h"^  =0(1)  (so  that   kll  B  ,11   is  bounded).  Theorem  2  is 
applicable,  provided  k  and  Eq  are  sufficiently  small.  One 
can  then  obtain  AP  sequence  solutions  and  carry  out  the 
convergence  argument  of  Section  10  (see  (10.17)  f f ) , 
with  the  only  difficulty  being  the  verification  that 

(11.30)  II  (B^-B^)yll^  =  O(h^) 

and  this  will  follow  immediately  if  one  assumes  a  chain  as 
in  Assumption  I  for  the  operator  L   introduced  in  (10.16). 
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Appendix.   Proof  that  (L^,L^)  Admissibility  Implies 
an  Exponential  Dichotomy. 

Here  we  would  like  to  give  a  simplified  proof  of  the 
proposition  that  (A, A)  admissibility  for  the  inhomogeneous 
equation 

(A.l)  y  =  Ay  +  f 

implies  an  exponential  dichotomy  for  the  homogeneous  equation 

(A.  2)  y  =  Ay 

It  is  assumed  of  course  that  A(t)  is  AP  although  we  will 
actually  do  the  L^  case.  The  proof  is  valid  only  if  the 
underlying  space  E  is  finite  dimensional  and  we  can  then 
regard  A  as  a  matrix  and  y  and  f  as  m  vectors,  where  m  is 
the  dimension  of  E.  The  reader  is  referred  to  Section  5 
where  an  entirely  similar  proof  for  the  difference  equation 
case  is  given. 

The  general  proof  is  given  in  M  &  S  (p.  344,  Theorem 
103. A) ,  but  this  proof  relies  strongly  on  preceding  material. 
A  finite  dimensional  proof  is  given  in  Coppel  (p.  134  ff) 
for  an  equation  on  [0,°°).   His  proof  can  be  extended  to  an 
equation  on  the  whole  real  axis,  but  the  proof  given  here 
is  simpler  and  more  in  keeping  with  the  ideas  of  Massera 
and  Schaffer. 

It  is  shown  in  M&S  (Theorem  10 3. A)  that  (A, A)  admissi- 
bility implies  (L  ,L  )  admissibility  and  this  will  be  our 
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starting  point.      We    thus    assume    that    for   every    f    in   L       ,     (A.l) 
has    a   unique   solution  y    in   L^   and   the   estimate 

(A.  3)  llyll       <   Kll  fll 

will  hold  for  a  certain  constant  K. 

Now  let  S   be  the  subspace  of  initial  data  which  gives 
rise  to  solutions  to  (A. 2)  that  are  bounded  for  t  s  [0,°°). 
Let  Sy      be  the  analogous  subspace  giving  rise  to  solutions 
that  are  bounded  for   t  G  (-°°,0].   We  must  have  S   n  s  =  {0} 
because  any   nontrivial,  bounded  solution  to  (A. 2)  would 
violate  the  uniqueness  requirement  of  our  definition  of 
admissibility.   Let  S _  be  any  subspace  complementary  to 
S,  ©  S2  so  that 

(A. 4)  E  =  S^  ©  S2  ®  S^ 

Let  PwPp,P^  be  the  associated  projections,  i.e.  P,  is  the 
projection  onto  S,  along  S„  ©  S^  ,  etc.   We  have 

(A. 5)  I  =  P^  +  P2  +  P3 

Pj^P.  =  0  ,   i  =  1,2,3;   j  =  1,2,3;   i  7^  j . 

We  will  first  show  S^  =  (O),  i.e.  P^  =  0.  This  will 
follow  from  the  variation  of  constants  formula,  just  as  in 
the  difference  equation  case.  If  f  has  compact  support  and 
y(t)  is  the  unique  L^  solution  to  (A.l),  then  since  y  is  a 
bounded  solution  to  (A. 2)  for  large  |t|,  we  must  have 
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(A. 6)   (a)    P^Yq   =   -    I  P2Y  """(s)  f(s)  ds 

0 


(b)  P^Yq    =  -  J  P3Y"^(s)  f(s)  ds 

0 
0 

(c)  P^Yq    =       f   P;l^~"^(s)  f(s)  ds 


0 
(d)    P^Yq   =       f   P3Y"-^(s)  f(s)  ds 


Here  Y(t)  is  the  fundamental  solution  to  (A. 2). 

Now  for  any  vector   z  e  E,  if  we  set  f  =  Xrr,  iiY(t)z 
then  (A.6b,d)  yield  immediately  P^z  =  0  and  we  can  conclude 

P^  =  0,  i.e. 

E  =  S   ©  S„ 
(A. 7)  ^  ^ 

I  =  P^  +  P2 

Using  (A.6a,c)  we  also  see  that  for  any  f  with  compact 
support  we  have 


(A. 8)  y(t)  =     G(t,s)  f(s)  ds 


where 

Y(t)P  ¥""""(3)  ,       t  ^  s 
(A. 9)  G(t,s)  =         _, 

-Y(t)P2Y   (s)  ,       s  ^  t 

We  will  next  use  (A. 3)  to  show  that  for  every  fixed  t, 
G(t,s)  is  in  L,  (-°°,°°)   as  a  function  of  s,  and  in  fact 
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(A. 10)  ds  llG(t,s)ll  <_   K 

—  oo 

where  K  is  given  in  (A. 3) .  This  is  a  simple  consequence  of 
using  (A. 8)  as  an  operator  on  L  (I),  where  I  is  any  compact 
interval  to  conclude  that 


(A.  11)  ds  llG(t,s)ll  £  K 

I 
and  then  letting  I  -^  (-=°/°°)  . 

We  now  will  show  that  (A. 10)  implies 
(A. 12)  llG(t,s)ll  =  0(1)  . 


To  see  this  we  simply  observe  that  Y   (s)  satisfies  the 


adjoint  equation 


(A. 13)  z'  =  -  Z  A(s) 


where  we  have  used  "'"  to  stand  for  d/ds.   Hence  fixing  t 
and  considering  the  region  s  <  t,  we  have 


(A. 14)  G(t,s)'  =  -  G(t,s)  A(s) 


Now  since  A  is  bounded  it  follows  that  G'  {t,s)  is  in  L,  ,  as 

a  function  of  s ,  on  the  interval  (-°°,t]  with  an  L,  norm 

which  can  be  bounded  uniformly  in  t.   This  of  course  means 

that  G{t,s)  must  approach  a  limit  as  s  ^  -°°  and  since  G(t,s) 

is  in  L   (as  a  function  of  s) ,  this  limit  can  only  be  zero. 

We  thus  have 

s 


(A. 15)        G(t,s)  =     dr  G' (t,r)  -  0(1)  ,   s  <  t. 


—  00 
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and  a  similar  argument  establishes  (A.  12)  for  s  >^  t. 

We  must  now  show  that  (A. 12)  can  be  replaced  by  an 
exponentially  decaying  factor,  and  to  do  this  we  consider 
perturbations  of  (A. 1)  and  (A. 2);  namely  for  small  e  >  0 
consider  the  equations 

(A. 16)  y  =  [A-eI]y  +  f 

(A. 17)  y  =  [A-£l]y 

If  e  is  sufficiently  small  we  will  have  (L  ,L  )  admissible 

-*  OO      QO 

for  (A. 16).   This  is  a  simple  consequence  of  the  contracting 
mapping  principle  analogous  to  the  argument  of  Section  4C  (see 
(4.9)  ff).  The  contracting  mapping  principle  will  also  show  that 
the  admissibility  bound  for  (A. 16)  will  be  given  by  (see  A. 3) 

(A. 18)  K^  =  — ^^-^-  IK,  Ul  1  ^n  ' 

^    (1-eK)  ^ 

if  e^    is  sufficiently  small.   Thus  the  admissibility  bound 
for  (A. 16)  can  be  taken  independent  of  e,  if  e  is  small 
enough . 

Now  the  fundamental  solution  to  (A. 18)  is 


(A. 19)  Y^  (t)  =  e  ^''^  Y(t) 


and 


if   Pw  P^  denote  the  projections  replacing  P,  and  P 


while  G^(t,s)  is  the  new  Green's  function,  then  the  analysis 
which  led  to  (A. 12),  and  which  depended  only  on  admissibility, 
is  equally  valid  here  and  we  can  conclude 
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(A.  20)  llG^(t,s)ll  =  0(1)  ;  "  G^  (t,  •  )  " -^  =  0(1)  , 

where  the  bound  in  (A. 20)  can  be  taken  independent  of  e 
(because  of  (A.18)). 

The  estimate  (A. 20)  means  in  particular  that  for  small 
e  there  is  a  constant  K_  independent  of  e  such  that  if  s>^  t, 

(A. 21)        llY(t)P2Y"-'-(s)ll  <K2e~^^^"'^^ 

and  it  is  only  necessary  to  show  that 

P^   =  P 
1      1 

(A. 22)  ; 

P    =  P 
^2      2 

To  prove  (A. 22)  we  let  S^  and  S2  replace  S,  and  S2. 

It  then  follows  from  (A. 19)  and  e  >  0 

•=;   C  c;^ 

1  -  1 
(A. 23)  ^ 

S2  c  S2 

and  (A. 22)  will  be  established  if  we  can  show  that  the  ranks 
of  the  stable  and  unstable  projections  are  unchanged  for 
small  e,  and  this,  in  turn,  will  follow  from 

(A. 24)  IIP^  -  P^ll  =  0(e)  . 

To  obtain  (A. 22)  we  set  H^(t,s)  =  G^(t,s)  -  G(t,s). 
Now  for  any  f  with  compact  support  we  let  y  and  y  be  the 
unique  bounded  solutions  to  (A.l)  and  (A. 16)  respectively. 
Ifw   =v-v  we  see  that  w   is  the  unique  bounded 
solution  to 
(A. 25)  w^  =  Aw^  -e  y^  , 
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and  we  can  conclude  from  (A. 3)  together  with  its  analogue 

for    (A. 16)    that 

(A.  26)  II  w    II       =   0(e)  II  fll 

'  £;     oo  oo 

Now  w      can   also  be  written    as 

oo 

(A. 27)  "e^t)    =  ds    H^(t,s)    d(s) 

—  OD 

and   using   the   same   argument   that    led   to    (A. 10)    we   see    that 

CO 

(A. 28)  I    ds    llH^{t,s)ll    =    0(e)     . 

—  oo 

(Note    that    the   bound   in    (A. 28)    will   be    independent   of    t, 
although  we  will   only    use    it    for    t  =    0.) 

Now   if   s   £  t  we    see    that   as    a   function   of   s,    H    (t,s) 
satisfies    the   equation 

(A. 29)  H^      =    -    H^A   +    e    G^(t,s)     . 

This   equation,    together  with    (A. 20)    and    (A. 28),    shows    that 

(A. 30)  I    ds    llH^     (t,s)ll    =    0(e) 

—  oo 

and  since  H  (t,s)  is  in  L,  as  a  function  of  s  we  can  write  for 

'  -  ^  t       ^ 

(A. 31)       H^(t,s)  =   I  ds  H^  (t,s)  =  0(e) 

—  CO 

and  if  we  set  t  =  0 ,  s  =  0_  ,  (A. 31)  yields  (A. 24),  which 
proves  one  part  of  the  requirement  of  an  exponential  dichotomy 
and  the  other  part  follows  on  replacing  -e  by  +e. 
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